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ABSTRACT 

We present new far infrared maps of the Small Magellanic Cloud (SMC) at 24, 70, and 160 /im 
obtained as part of the Spitzer Survey of the Small Magellanic Cloud (S'^MC. feolatto et al.ll200^ . 
These maps cover most of the active star formation in the SMC Bar and the more quiescent Wing. 
We combine our maps with literature data to derive the surface density across the SMC. We find 
a total dust mass of Mdust =3 x 10^ M0, implying a dust-to- hydrogen ratio over the region studied 
of logiQ D/H = —2.86, or l-to-700, which includes H2. Assuming the dust to trace the total gas 
column, we derive H2 surface densities across the SMC. We find a total H2 mass Mh2 = 3.2 x lO"^ M© 
in a distribution similar to that of the CO, but more extended. We compare profiles of CO and 
H2 around six molecular peaks and find that on average H2 is more extended than CO by a factor 
of ^ 1.3. The implied CO-to-H2 conversion factor over the whole SMC is Xqo = 13 ± 1 x 10^^ 
cni^^ (K km s^^)^^. Over the volume occupied by CO we find a lower conversion factor, Xqo = 
6 ± 1 X 10^^ cm^^ (K km s~^)^^, which is still a few times larger than that found using virial mass 
methods. The molecular peaks have H2 surface densities Eh2 ~ 180 ± 30 Mq pc~^, similar to 
those in Milky Way GMCs, and correspondingly low extinctions. Ay ~ 1 - 2 mag. To reconcile 
these measurements with predictions by the theory of photoionization-regulated star formation, which 
requires Ay ~ 6, the GMCs must be ~ 3 times smaller than our 46 pc resolution element. We find 
that for a given hydrostatic gas pressure, the SMC has a 2 - 3 times lower ratio of molecular to atomic 
gas than spiral galaxies. Combined with the lower mean densities in the SMC this may explain why 
this galaxy has only 10% of its gas in the molecular phase. 

Subject headings: ISM: dust, ISM: molecules, Galaxies:Magellanic Clouds, Galaxies:ISM, Infrared: 
ISM 



1. INTRODUCTION 

Dust grains play a central role in the formation of stars 
from atomic gas. Dust shields molecular gas from dissoci- 
ating radiation and allows it to cool and condense to the 
densities necessary to form stars. The surfaces of dust 
grains serve as the sites of molecular hydrogen forma- 
tion. In low metallicity galaxies, an underabundance of 
dust may result in more intense radiation fields, different 
thermal balance in the interstellar medium (ISM), and 
different structure for giant molecular clouds (GMCs). 
However, there are only a few low metallicity systems 
in which atomic gas, molecular gas, dust, and star for- 
mation tracers have all been mapped at resolution suffi- 
cient to resolve individual star-forming regions or GMCs 
(< 100 pc). 

In this paper we present new far infrared (FIR) images 
of the SmaU Magellanic Cloud (SMC) obtained using the 
Spitzer Space Telescope as part of the S pitzer Survey o f 
the SmaU Magellanic Cloud fS^M CBolatto et al.ll20'06^ . 
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We combine these maps with previous FIR, H I and CO 
data to construct a complete picture of the star-forming 
ISM in this nearby low metallicity galaxy. We focus on 
three questions: 1) what is the mass of dust in the SMC?, 
2) how much molecular hydrogen (H2) is in the SMC? 
and 3) how does the distribution of H2 in the SMC — 
particularly its structure and relation to CO emission — 
compare to that in the Milky Way? In this introduction 
we describe the SMC and expand on these questions. 

The SMC is the nearest actively star-forming, low 
metallicity {Z < 1/10 Zq) galaxy. It is therefore good 
location to study the relationship between gas, dust, 
and star formation at low metallicity. It forms stars 
at ~ 0.05 Mq yr-i fWilkc et al. 20041) — a normal- 
ized rate that is comparable to that of the Milky Way 
(the dynamical mass of the SMC is 2 x 10^ Mq — 
IStanimirovic. Stavelev-Smith. fc Jones! l2004l ) . It is gas 
rich, with ft! 4 x 10^ M© of atomic gas in an extended 
distribution (jStanimirovic et al. 1999) and a comparable 
mass of stars, 3 x 10® Mq (estimated from the DIRBE 
2.2 /im data). The SMC is also unenriched in heavy el- 
ements, with an oxygen a bundance just above a tenth 
solar (12-H og(0/H) = 8 .0,'Dufour"l984'). At a distance 
of 61.1 kpc ()Hilditch et al. 2005; Keller & Wood 2006), 
even modest telescopes can observe the SMC with ex- 
cellent spatial resolution, and at a Galactic latitude of 
b = —44° it is obscured by little foreground material. 

Several recent studies have used FIR emission to study 
the dust in the S MC, finding eviden c e for a low dust-to- 
gas ratio, DGR. IStanimirovic et aTl (|2000( ) modeled HI- 
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RAS data (a high resolution IRAS product) and found 
that the SMC had an unexpectedly low dust mass of 
only 1.8 X 10^ Mq. The corresponding PGR is two or- 
ders o f magnitude below the Galactic value. iBot et ahl 
12004') studied the diffuse gas using ISO and IRAS data, 
and found that the emissivity of diffuse SMC gas is 
30 times lower tha n that of high latitude Galactic gas. 
iWilke et all ((200l studied emission integrated over the 
whole SMC, also using ISO data. They fit the dust using 
several black bodies and found a much higher total dust 
content than previous studies, 7.8 x 10^ M©. Using this 
result, they derived a DGR only a few times lower than 
the Galactic value. 

Only 1.7 X 10^ M© of the dust mass obtained by 
IWilke et all ()2004f ) would be estimated from their FIR 
data alone. The rest is in a cold component that emits 
at millimeter wavelengths , with a temperature of ~ 10 K. 
IStanimirovic et ahl (|2000f ) also recognized the possibihty 
of a large mass of cold dust from the 140 and 240 /^m 
DIRBE data, which show higher intensities than those 
allowed by a single popula tion of Galact i c dust that fits 
the IRAS measurements. lAguirre et al.l ()2003D demon- 
strated that this excess extends out to 1.2 mm by mea- 
suring the spectrum of the SMC using the TopHat bal- 
loon. They also showed, however, that the measurements 
can be fit by a single dust population if /3 = 1 for the 
dust emissivity law, r oc A^*^. 

In addition to providing information about the dust, 
the FIR emission can be used to trace molecular gas. 
Clouds of molecular hydrogen are the sites of star for- 
mation, but II2 cannot be observed directly under the 
conditions found in GMCs because it lacks a permanent 
dipole moment and T is too low to excite the quadrupole 
transitions. Observations of tracer molecules, particu- 
larly CO, often substitute for direct observations of H2. 
In low metallicity galaxies, the lack of heavy elements 
and the lower DGR may alter the structure of molec- 
ular clouds and affect the relationship between tracer 
molecules and H2, leaving the amount of molecular gas 
in these systems uncertain. Unlike the CO-to-Il2 ratio, 
the dust-to- hydrogen ratio, -D/iJ, can be measured away 
from molecular clouds using H I and FIR emission. Un- 
der the assumption that D/H varies slowly, FIR emis- 
sion may provide a better estimate of the amount of II2 
in these systems than CO observations. 

The molecular hydrogen content of the SMC may be 
even l e ss wel l known than its dust content. iMizuno et al.l 
(|2001l I2006t l surveyed the CO J= 1 ^ line across 
most of the SMC using the NANTEN telescope and 
found a luminosity of 1 x 10^ K km s""'^ pc^. They 
used the virial mass method to derive a CO-to-Il2 con- 
version factor, Xcoi of 1-5 xlO^^ cm~^ (K km s~^)~^, 
about 10 times the Galactic value (A^co*^'^' « 2 x 
10^° cm^^ (K km s^^)^^). This implies a t otal molec- 
ular mass of A/moi ~ 2 - 10 x lO*' M©. iBlitz et al.l 
([2006) corrected these data for resolution and sensitivity 
effects and found a value of Aco towards the low end of 
this range, 1 - 1.5 x 10^^ cm~^ (K km s~^)~^, implying 
-^moi ~ 2-3 X 10^ M©. Consi dering individual re gions 
with better spatial resolution, ' Rubio et al.l (|1993f ) and 
[Bplatto ct al. (2003) both found conversion factors a few 
times the Galactic value, though dependent on the size 
of the structure studied. 



Modeling of the FIR can be used to derive the mass 
and structure of H2 without relying on the assumptions 
that CO and II2 have the same structure and that GMCs 
are in virial equilibrium. If the gas and dust are well 
mixed and properties do not vary too dramatically be- 
tween the molecular and atomic ISM, the dust acts as 
an optic a.lly thin trac e r of th e total, H I-I-H2, gas surface 
density. I Dame et al.l (|200lD demonstrated this method 
to work comparing IR A S and CO in the Milky Way (see 
also lReach et al.lll998D . Ilsraell (|1997D used IRAS to study 
the SMC and found a total H2 content of 7.5 x 10^ M© 
and the CO-to-H2 conversion factor to be 60 times the 
Galactic value. In this paper we use a similar method to 
construct II2 maps with 46 pc resolution and sensitivity 
to H2 surface densities as low as ^ 20 M© pc~^ (Sc). 

Using the II2 map, we test three predic tions regarding 
molecular gas in low metallicity systems. iMcKed i) 19891 ) 
argued that photoionization by far ultraviolet photons 
sets the level of ionization inside GMCs and thus regu- 
lates their structure. Magnetic fields couple to the ion- 
ized gas and support cores against collapse. Thus, the 
level of ionization is tied to the rate of star formation in 
the cloud. Low mass star formation provides feedback 
that keeps the mean extinction, and with it the ioniza- 
tion level, across the cloud _at a particular value. In the 
Milky Way this value is Ay ~ 4 - 8, independent of 
cloud inass. I n low metallicity systems, with low DGRs, 
iMcKed ([1981 predicts GMCs wiU require higher gas sur- 
face densities, Ehi+h2, in order to achieve the same mean 
extinctions found in the Milky Way. For the low DGR 
in the SMC, this impfies GMCs w ith Shi+H2 > 1000 
M© pc-2 (see also lPak et al.l 119981 ). Second, CO in the 
SMC may be underabundant or absent in regions where 
II2 survives because the II2 self-shields but this does not 
stop photodiss ociation of the CO (e.g..lMalonev Sz BlacS 
11988). Third, IBlitz k Rosolowskvl (j200j) showed that 
there is a good correlation between the midplane hy- 
drostatic pressure and the ratio of molecular to atomic 
gas along lines of sight towards spiral galaxies. We test 
whether this correlation extends to the low metallicity, 
low pressure environment found in the SMC, or whether 
metallicity acts as a "third parameter" in the relation- 
ship between pressure and molecular gas fraction. 

This paper is organized as follows: in Sj2] we discuss 
the new Spitzer maps and the literature data we use 
to supplement them. We describe our foreground sub- 
traction and discuss uncertainties in the data. In SjSjwe 
present the new maps and literature data on a common 
astrometric grid. We construct maps of dust and II2 
and present an updated FIR spectral energy distribution 
(SED). In [|4|we consider the relationship between dust 
and gas, including the dust-to-gas ratio, the CO-to-Il2 
conversion factor, the structure of molecular clouds, and 
pressure and extinction in the star-forming regions. In 
^we present our conclusions. 

2. DATA 

To meet the science goals described in the introduc- 
tion, we require FIR images of the SMC at several wave- 
lengths on a common astrometric grid with good relative 
calibration. These allow us to derive maps of the dust 
and II2 content in §|3|and add new points at 24, 70, and 
160 fim to the integrated SED of the SMC. We start with 
the new Spitzer images of the SMC from S^MC at 24, 70, 



S3MC: FIR in the SMC 



3 



and 160 /ini and supplement these data with images at 60 
and 100 /xm from IRAS and at 170 ^m from ISO. We use 
data from DIRBE from 25 - 240 /im as th e basis for the 
flux c ahbration and consider TopHat data (jAguirre et al.l 
l2003l ) when constructing the integrated SED. The wave- 
lengths, telescopes, limiting resolution, and references for 
the data used in this paper are listed in Table [T] In the 
rest of this section we describe how we isolate the emis- 
sion from the SMC in the FIR images and estimate the 
uncertainties in the final maps. 

2.1. Color Corrections and Flux Definitions 

Throughout this paper, we quote monochromatic flux 
densities assuming a. F^, (x spectral energy distribu- 
tion across the relevant bandpass. This is the conven- 
tion for IRAS, IRAC, and DIRBE, but differs from the 
MIPS standard, which is to assume a Rayleigh- Jeans tail 
{F^ oc z/^) across the bandpass. We apply color correc- 
tions to the MIPS data at 24, 70, and 160 /im by dividing 
by 0.961, 0.918, and 0.959, respectively, before any pro- 
cessing to bring them to the IRAS deflnition. For most 
of the work in this paper, this choice is arbitrary. 

2.1.1. Bright Pixel Removal 

Before other processing we remove pixels with unex- 
pectedly high intensities from the 24, 70, and 160 /im 
maps (> 100, > 200, and > 250 MJy ster~^, respec- 
tively). At 24 /im, these pixels correspond to point 
sources and are well over the saturation limit, at 160 
/xm these are simply artifacts, and at 70 //m they are a 
mix of the two. Removing these point sources lowers the 
integrated flux in the foreground subtracted 24 fim map 
by « 15%. The effect on the flux at 70 /im is negligible 
(~ 1%). We interpolate to fill these pixels using nearby 
data when placing the maps at a common resolution by 
taking the beam-weighted mean of all nearby unflagged 
data. 

2.2. Foreground Removal 

The intensity in the FIR maps is the sum of emission 
from the source and several foregrounds, 

^FIR = ^Zod + -^Gal+CIB + ^Source • (1) 

Izod is the intensity of zodiacal light, which dominates 
the foreground at shorter wavelengths (12 - 70 /im). 
Icai+ciB is Galactic dust emission and a small diffuse 
cosmic infrared background term. This is the primary 
foreground at longer wavelengths (100 - 240 /im). 

2.2.1. Zodiacal Light 

The HIRAS, IRIS, ISO, and DIRBE data have all 
been corrected for zodiacal light. We estimate the con- 
tribution of zodiacal light to th e MIPS maps using the 
model of interplanetary dust by iKelsall et al.l (|l998r ) im- 
plemented in the Spitzer Observation Planning Tool. A 
planar fit derived to this model near the SMC for our 
observations is: 

Izod = lo.zod [0.006(/ - IsMc) - 0.003(6 - bsMc) + 1] 

(2) 

where Iq zod is the intensity of zodiacal light at the nom- 
inal center of the SMC {I = 302°. 7969, b = -44°.2992) 



on the date of the observations. At 24, 70, and 160 /im, 
respectively, /o,zod is 20.1, 5.3, and 0.9 MJy sterol (98%, 
36%, and 6% of the total flux, respectively). The varia- 
tion in Izod over the several days during which the ob- 
servations were taken is minimal. 

2.2.2. Galactic Foreground and Cosmic Infrared Background 

At 100 /im and longer wavelengths the Galactic ISM 
contributes the largest foreground. We apply a similar 
treatment to lBot et al.l |2004), and estimate the Galactic 
foreground along lines of sight towards the SMC using 
the Galactic H I column density, which is clearly sepa- 
rated from the SMC H I column density by its velocity. 
We use the Galactic H I column density o bserved at ^ 
0.5° re solution from the southern survey bv lBajaja et al.l 
(|2005f ). For the IRIS 100 /im map, we us e the fit of 
FIR in tensity to H I column density by .Boulanger et al.l 
dUll). At 160 /im and 170 /im, we interpolate from the 
fits at 140 and 240 /im by iBoulanger et all (|1996l ) into 
these bands using their best fit temperature T = 17.5 K 
and emissivity /3 = 2. The median intensity of the Galac- 
tic foreground plus cosmic infrared background (GIB) is 
1.8, 4.7, and 4.8 MJy ster"! at 100, 160, and 170 /im, 
respectively. 

For the 60 and 70 /im m aps, the foregrou i id du e to 
Galactic dust is small, but IBoulanger et al.l ()1996t ) do 
not provide fits for these bands. We use estimates of 
the emissio n of local dust per unit H I column density at 
60 /im from lDesert et all ()1990. ') and interpolate to 70 /im 
using the average colors in the map. This calculation is 
approximate, but the median intensity of this foreground 
at these wavelengths is small: 0.5 MJy ster^^ and 1 MJy 
ster~^, or < 20% of the zodiacal light foreground. 

We do not subtract any foregr ounds from the HIRAS 
maps; IStanimirovic et al] ()2000[ ) define regions of sky 
near the SMC as blank and subtract a single foreground 
that includes zodiacal light, the Galactic foreground, and 
the GIB. Because we do not have an independent esti- 
mate of the zodiacal light for these data, we cannot break 
this foreground into its individual components. Instead, 
we bring the HIRAS data onto a consistent foreground 
subtraction by comparing them to the foreground sub- 
tracted IRIS data at a common resolution. We find that 
offsets of 0.75 MJy ster"! and 1.5 MJy ster^^ at 60 and 
100 /im, respectively, are needed to bring the HIRAS zero 
level into agreement with IRIS. After this correction, the 
agreement between the two IRAS data sets is excellent. 

2.2.3. Foreground Subtracting the DIRBE Data 

The DIRBE Zodiacal-Subtracted Mission Average 
(ZSMA) maps appear to have a small amount of resid- 
ual zodiacal light, ~ 1 MJy ster~^ judging from a plot 
of intensity as a function of H I column over the whole 
sky. Therefore, we treat the DIRBE foreground subtrac- 
tion separately. We derive DIRBE intensities as a func- 
tion of H I column density near the SMC (within 10° of 
the optical center) but offset enough that we see only 
Galactic FIR emission (more than 3° away from the op- 
tical center). Then, for lines of sight towards the SMC 
we measure the Milky Way H I column density as above 
and subtract the average intensity of lines of sight with 
the same H I column density away from the SMC. This 
is our estimate of the DIRBE foreground, which should 
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TABLE 1 
Data Sets Used in This Paper, 



Nominal Wavelength 


Data Set 


Limiting 


Spatial 


Reference 






Resolution 


Resolution^ 





24, 70, 160 

25, 60, 100, 140, 240*^ 
25, 60, 100 

25, 60, 100 



MIPS (Spitzer) 
DIRBE (CORE) 
HIRAS (IRAS) 
IRIS (IRAS) 
ISOPHOT (ISO) 
ATCA + Parkes 

NANTEN ^ 

" At our adopted distance of 61.1 kpc jHilditch et al . 2005: Keller & V \^odf 
[2006!). The DIRBE zodiacal light subtracted mission average (ZSMA) data 
are the basis of our photometric calibration in this paper. 



170 
HI 

CO 



40" 
0°.7 
1'.7 
4'. 3 
1'.5 
1'.5 
2'. 6 



12 pc 
750 pc 
30 pc 
76 pc 
27 pc 
27 pc 
46 pc 



Bolatto et al. (2006) 
H auser et al. (1998) 
Stanimirovic et al. (2000) 
Miville-Dcscliones & Lagache (200^ 
Bot et al. (2004) 
Stanimirovic e t al. Q999) 
Mizuno et al. (2001. 2006') 



account for Galactic emission, the GIB, and residual zo- 
diacal light. 

The DIRBE 25 /im data is overwhelmingly contami- 
nated by zodiacal light. We fit a plane to the nearby 
data and subtract this foreground before using it. 

2.2.4. The MIPS 24 ^im Map Foreground 

In the case of the 24 /xm map there are large regions 
of the map free of diffuse emission. Because the zodia- 
cal light foreground is large and potentially uncertain at 
24 /im and we do not have an estimate of the Galactic 
foreground, we use blank regions in the Spitzer map to 
estimate the residual foreground. We measure an addi- 
tional foreground with a median value of 2.5 MJy ster~^ 
above and beyond the original zodiacal light foreground 
(i j2.2.ip . This value is probably too large to be due to 
the to Galactic foreground alone and may represent a 
failure in the zodiacal light model or a calibration error 
for MIPS. The presence of blank sky in the 24 /xm map 
makes it easy to remove, regardless. 

2.2.5. Consistency and Calibration 

We check the calibration of the high resolution data 
against DIRBE and find generally good agreement. The 
fluxes from each map integrated over the whole SMG 
agree (see Table O. The exception is the MIPS 160 
/xm map, which contains substantially more flux than 
we estimate based on the DIRBE data. A central goal 
of the DIRBE missio n was to obtained a bsolute bright- 
ness maps of the sky (iHauser et al.lflQQSl) , and that data 
represents a more reliable measurement of the flux than 
our MIPS map. Therefore, we calibrate the 160 /xm to 
match DIRBE by rescaling the map. Appendix |^ de- 
scribes how we carry out this calibration in detail. The 
MIPS 70 /.tm map is consistent with the data at 60 yum 
and 100 /xm, but because the SED has a very steep slope 
over the MIPS 70 /im bandpass, this consistency is not a 
particularly strong constraint. It is not possible to per- 
form the same sort of detailed calibration we adopt at 
160 /xm for the 70 /xm map. 

2.2.6. Saturation 

Saturation may represent a concern over part of our 
160 /Ltm map. The bright star-forming region in the 
southwest of the galaxy contains a significant region with 
brightness greater than 50 MJy ster~^, the nominal sat- 
uration limit for our medium speed scan maps. A few 
other star-forming regions (N 76 to the north and N 83 
in the east) also contain small regions with intensities 



above 50 MJy ster~^ In total, 13% of the flux and 2.5% 
of the area are in regions with intensities above the sat- 
uration limit. In these regions, the response of detector 
becomes nonlinear and we have an additional gain un- 
certainty. The figure in Appendix A shows that in the 
southwest of the bar — the source of the lines of sight 
with intensities > 25 MJy ster~^ at the DIRBE resolu- 
tion — the MIPS data are ~ 20 - 30% lower than the 
DIRBE data. Outside this region, the MIPS and DIRBE 
data are well matched after correction. 

2.3. Uncertainties in the Final Maps 

We estimate the statistical uncertainties for the FIR 
data from the scatter about the median in regions of low 
intensity emission. We fit a Gaussian to the residuals 
about the local median taken over at 4' x 4' region, so 
that for each point Iresid = Imap - Imedian (±2') with 
Imedian (i2') the median over a 4' x 4' box centered on 
that point. We conduct these measurements in regions 
of low intensity emission outside the main body of the 
galaxy to minimize the effect of structure on the mea- 
surement. However, emission with small scale structure 
and artifacts present in the MIPS and ISO maps con- 
tribute to the noise values. These are legitimate sources 
of uncertainty, because they affect measurements. 

Table [2] gives the results of our noise estimates, and 
Figure [1] shows a map and histogram of deviations about 
the median for the MIPS 160 /im data. Two artifacts of 
small scale structure are visible: the tail of low, but non- 
zero, values towards negative intensities and the slight 
shift of the center of the Gaussian towards negative in- 
tensities. Table [2] gives the noise for each map at the 
native and IRIS resolution. For most maps, the noise 
scales approximately as expected with averaging. Devia- 
tions from this scaling may be the result of artifacts such 
as the visible striping at low intensity in the MIPS data. 

Because we calibrate our data against DIRBE, the ab- 
solute uncert a inties in DIRBE are relevant to our data. 
IHauser et al.l (|1998D gives the absolute gain uncertainty 
for DIRBE at 25, 60, 100, 140, and 240 /im as 15.1, 10.4, 
13.5, 10.6, and 11.6 %, respectively. The uncertainties 
in the absolute offset of the DIRBE data are all small 
enough to be negligible for our purposes. As a result of 
gain uncertainties for DIRBE, the minimum systematic 
uncertainty in any absolute flux or intensity in this paper 
is ~ 10%. 

The systematic uncertainty associated with the fore- 
ground subtraction is ±1 MJy ster"^. This is the ap- 
proximate magnitude of the scatter about our best fits 
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TABLE 2 
Random Uncertainties in FIR Maps 



Data Set 








Native Resolution^ 


IRIS Resolution'' 




MJy ster~^ 


MJy ster~l 


24 /im'» MIPS (Spitzer) 


0.02 


0.01 


60 /im HIRAS (IRAS) 


0.07 


0.02 


60 /xm IRIS (IRAS) 


0.02 


0.02 


70 Atm'' MIPS (Spitzer) 


0.34 


0.06 


100 Atm HIRAS (IRAS) 


0.10 


0.05 


100 Atm IRIS (IRAS) 


0.05 


0.05 


160 Atm MIPS (Spitzer) 


0.58 


0.10 


170 Atm ISOPHOT (ISO) 


0.22 


0.05 



24 and 70 nm data convolved to 160 fim resolution.'' Uncer- 
tainty after convolution to the IRIS (4') resolution. 




Fig. 1. — (top) The map used to measure the noise at 160 A*™. 
The map shows fluctuations of the 160 Aim MIPS map about the 
local median (over a 4' X 4' box centered on the pixel). Regions 
with bright emission or no data are blanked out and appear black. 
Artifacts, particularly striping along the direction of the map, con- 
tribute to the noise, but very little extended structure is visible. 
The extent of the MIPS 160 A^m map is shown by the white bor- 
der, (below) The histogram of the residuals about the median. The 
dashed lines shows the best fit Gaussian, which has crieo = 0.58 
MJy ster-l. 

of DIRBE intensity to H I column density (~ 0.2 MJy 
ster-i at 100 ^m, ~ 1.2 MJy stei~^ at 240 ^m); it is 
also about the magnitude of the offset needed to bring 
IRIS and HIRAS into rough agreement (~ 0.75 and 1.5 
MJy sterol). 

3. RESULTS 

In this section we present maps of FIR emission and 
neutral gas from the SMC. We use these maps to de- 
rive the dust and molecular gas surface density and we 
present an updated SED of the SMC that includes the 
new 24, 70, and 160 fim data. 

3.1. FIR Maps of the SMC 



Figure [2] shows the relative extent of 160 iim (red), 24 
Atm (green) , and 8 /x m (blue) emission in the S^MC maps 
(|Bolatto et al.ll2006[ ). The 160 ^m emission is bright over 
the whole SMC Bar in the west and in the N 83/N 84 
star- forming region far to the east. This emission comes 
mostly from big grains (a « 15 - 110 nm) with T ~ 20 K 
with a small (few percent) contribution from the [CII] 158 
/im line (see Appendix E]) . Every region in the S MC 
map shows some 160 /im emission including the SMC 
Wing, which is roughly perpendicular to and east of the 
bar and contains more diffuse gas. By contrast, the 24 
/zm emission and the 8 /xm emission are both confined to 
regions of active star formation and point sources. Much 
of the 24 /im emission comes from stochastically heated 
very small grains, though very hot ('-^ 100 K) large grains 
may also make a contribution. In environments similar 
to the Solar Neighborhood, most strong 8/im emission 
appears to originate from polycyclic aromatic hydrocar- 
bons. However, very small grains may produce emission 
in this band even if PAHs are absent. All three bands 
are brightest towards the N 66 star-forming region in 
the southwest part of the galaxy. This is the brightest 
star-forming r egion in the SMC, home to 60 O stars 
(|Massev et al.l ll989'). 

Figures [3] and [4] show FIR emission and neutral gas in 
the SMC. FigureOshows the MIPS and IRIS images after 
foreground subtraction and processi ng. Figure [H shows 
H I co lumn density for the SMC from lStanimirovic et alj 
(| 1999( 1 an d the CO inte r isity map fr om the NANTEN 
survey by iMizuno et all ([2?)0ll [200l . The CO map is 
masked using a threshold contour of 0.3 K km s^^ (30"), 
which is then expanded by 1' in all directions. 

3.2. Dust Content of the SMC 

Figure [5] shows dust surface density across the SMC. 
We describe how we construct the map in Appendi x [Bl 
Briefly, we use the models by iDale fc Heloul (|200l ) to 
derive the dust mass surface density across the SMC. 
These models assume the dust along each line of sight 
to be heated by a power law distribution of radiation 
fields; we estimate the power law index from the 100-to- 
160 /zm color. Appendix |B] also gives the uncertainties 
associated with the dust map: a small < 10% statistical 
error, an additional 20% systematic uncertainty from our 
foreground subtraction, and a factor of 2 to 3 systematic 
uncertainty associated with the choice of emissivity and 
model are the most important. Over the region covered 
by the MIPS 160 nm map, we find a total dust mass of 
Mdust =3 X 105 M0. 

The SMC must have Mdust > Mq of dust; this 
value comes from assuming that a single population of 
grains is responsible for the emission in the IRAS 60 and 
10 liia bands. This is approximat ely the value obta ined 
bv IStanimirovic et all 12000) and ISchweri'^ (|1988D . In 
this case, the emissivity will need to have a wavelength 
dependence of /3 1 in order to match the emission 
observed at 140 /tm and beyond. 

Several estimates lead to an upper limit to the dust 
mass of Mdust ^ 10^ Mq. A Galactic dust-to-gas ratio, 
DGRmw ^ 10~^, apphed to the mass of H I in the SMC 
yields Mdust < 4 x 10^ Mq x 10"^ « 4 x 10^. It is well 
established that the DGR in the SMC is lower than the 
Galactic one, so this represents a firm upper limit. The 
heavy element abundance in the SMC is « 0.15 that of 
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Fig. 2. — A three color image showing the S^^MC maps l IBolatto et al.ll2003 ) of the SMC at 160 fim (red), 24 /xm (green), and 8 fj.m 
(blue). The 160 fim emission comes from big grains and extends across the whole SMC. Diffuse 24 fim emission, from hot very small grains 
(VSGs), and the 8 ^m emission, from PAHs, are both confined to the regions of active star formation. 



the Galaxy (measured from HII regions, IDufoud[l98l 
and the total H I conten t of the SMC is 4 x 10^ Mq 
(jStanimirovic et aLlll999f ). so the total mass of heavy el- 
ements in the ISM of the SMC is ~ 1.2 x 10^ M©. Be- 
cause the mass of dust cannot exceed the mass of heavy 
elements, this represents another, similar, upper limit to 
the dust content. 

The observed sub- mm/millimeter SED also constrains 
the mass of dust. We subtract the best fit {(3 — 2) 
single population model from the integrated spectrum 
and fit a single population to the residuals at long wave- 
lengths. Assuming this emission comes from cold dust 
with a Galactic emissivity (K250 — 8.5 cm^ g"-'^ and 
/3 = 2), we find a cold dust population with T « 7 K and 
Mdust ~ 2 X 10^ Mq. This is very close to the minimum 
equilibr ium teinperat ure possible for interstellar grains 
(~ 6 K, lPurcel]|[l969[ ) so this represe nts the largest dust 
mass per mitted by the SMC SED. IStanimirovic et aTl 
(j2000( ) and lWilke et all (|2004j ) both arrive a t simil ar val- 
ues for a possible, but not necessary (see ij3.3.ip . cold 
popi ilation, ^ 10^ M© and .7 x 10^ M©, respectively 
(see IStanimirovic et all l2000l . for the tradeoff between 
Tcoid and Mcoid)- 

3.3. The Integrated FIR SED of the SMC 



Table [3] gives the flux density integrated over the maps 
in the wavebands used in this pape r. We also give the 
range of literature values, taken from lWilke et al. (2004), 
who present an excellent compilation of the FIR SED of 
the SMC. Figure [6] shows log-log and log- linear plots of 
flux density as a function of wavelength. In both Table [3] 
and Figure [6l the area considered is the whole extent of 
the maps in Figure[3l not just the smaller area covered by 
the MIPS maps. Because this is an integrated SED over 
a large area, we are able to inclu de points from DIRBE 
and TopHat (jAguirre et al.ll2003f ). We scale the latter by 
0.9 to account for the difference between our area and the 
larger area studied by TopHat; this is the mean ratio of 
DIRBE fluxes over the TopHat area to the DIRBE fluxes 
over our own region. In deriving MIPS fluxes, we flU in 
the areas not covered by the MIPS maps with DIRBE 
data interpolated to the MIPS wavebands. These filled 
in regions account for ^ 20% of the total flux (the MIPS 
maps alone yield integrated fluxes of 310, 10000, and 
15000 Jy at 24, 70, and 160 /zm). 

Figure [6] and Table [3] show the photometric consis- 
tency discussed in ij2]above. The data for each waveband 
agree within the un certainties. Figure [5] also shows the 
iDale fc Heloul (2002) models (solid line) and best flt sin- 
gle populations with (3—2, 1.5, and 1 (dotted, dashed, 
and dash-dotted lines, respectively). All four models are 
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Fig. 3. — Maps of the SMC in six FIR bands arranged in order of increasing wavelength: (top left) the MIPS 24 fim map, (top right) 
the IRIS map, (middle left) the MIPS 70 fj,m map, (middle right) the IRIS 100 fim map, (bottom left) the MIPS 160 fim map, and 

(bottom right) the ISO 170 /im map (Bot et al. 2004). The MIPS maps are shown at resolutions of 40" with the extent of the map shown 
by a black contour. The IRIS maps have 4' resolution. The ISO map has 1'.5 resolution. All six maps have been corrected for zodiacal 
light and the Galactic foreground. The stretch is linear with the maximum for all maps selected to be 10 times the average intensity over 
a common area centered on the SW star-forming region. 
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Fig. 4. — The ISM in the SMC: (left) the H I at 1'.5 resolution l|Stanimirovic et al.|[l999l ^. and (right) CO intensity at 2'. 6 mapped by 
NANTEN lIMizuno et al...200L. .2009) . The black contour shows the extent of the NANTEN survey. 
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Fig. 5. — The dust mass surface density map for the SMC at 4' resolution. The black line shows the extent of the MIPS 160 fim map. 



derived from only the integrated 100 and 160 /im data 
to isolate emission from big grains. The best fit temper- 
ature for the big grains over the whole SMC is 20.6 K, 
somew hat higher than the T = 17.5 K found for Galactic 
cirrus (jBoulanger et al.l 11996 ^. The best fit value for a 
over the whole SMC is 2.4, implying cirrus-like medium 
with a mean radiation field close to the Solar Neighbor- 
hood value. 

3.3.1. The Long Wavelength Emission 

The SMC emits more at submillimeter and millime- 
ter wavelengths than one would estimate from the 
FIR emission assuming Galactic dust ((3 = 2). This 
may t ypical for a dwarf irregular galaxy ( Bolatto et al.l 
[2OOOI: iLisenfeld et an[200l IGalhano et al.l 120051) . This 
long wavelength emission may indicate the presence of 
cold dust invi sible in the FIR (Stani mirovic et al]l2000t 
lAguirre et cd] [2003; Wilkc et al.. ,2004 ): therefore, its in- 
terpretation is relevant to this work. If the long wave- 



length emission is due only to cold dust with a Galactic 
/? = 2, then it could have Tcoid ~ 7 K and Mcoid ^ 
2 X 10^ Mq, representing the dominant component of 
both d ust and heavy elements in the SMC. [Galliano et al.l 
(|2005l ) argued for this in several galaxies similar to the 
SMC: that very cold (5 to 9 K) dust makes up 40 to 80% 
of the dust mass. 

However, Figure [6] shows that the long wavelength 
emission can also be explained by an emissivity power 
law index shallower than (3 = 2 and little or no cold 
dust. The best fit to the FIR adopting /3 = 1.5 yields 
fluxes within 2cr of every TopHat data point, though still 
systematically low. The best fit /3 fr om 100 iim to 1 mm 
with a single population is P ~ I (jAguirre et al.ll2003l . 
find f3 = 0.91). Values of /3 = 1 and 1.5 are within the 
range of plausible astrophysical values, being appropriate 
for amorphous carbon and silicate grains, respectively. 
Because the dust in the SMC is believed to be dominated 
by silicate grains, a value of /? = 1.5 may be more appro- 
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priate than (3 = 2.0. Indeed, th e best fit (3 to the long 
wavelength values Kgbs found by IWeingartner &: Draind 
(l200ll):lLiirDraind (l200l L[200l for SM C grains beyond 
- SOOAtm is 1.65. Further. [Alton et all (|2004f) argue that 
/? = 1.5 describes emissivity determinations from the lit- 
e rature as well or better t han (3 = 2.0. 

iLisenfeld et"all (|2002l |2005[ ) observed similar long 
wavelength emission in NGC 1569, an irregular galaxy 
similar to the SMC in mass. They suggest that the 
small, hot grains responsible for the 24 - 70 /zm emis- 
sion may also cause the lon g wavelength emiss ion. Us- 
ing a modified version of the lDesert et al.l (|1990( ) models, 
they showed that an SED dominated by very small grain 
(VSG) emission (with no PAHs) can fit the data as well 
as a model containing cold big grains. They argue that 
NGC 1569, a metal-poor galaxy with a strong ISRF, is 
an unlikely locale to find a large amount of dust at un- 
usually low temperatures. While such dust might hide 
in a few molecular regions across the galaxy (we identify 
likely spots for this in the SMC below) there is a dearth 
of places to hide the truly large amount of dust needed 
to account for the long wavelength emission. Further, if 
such a reservoir of hidden does does exist, it necessarily 
means that the metallicity of these galaxies has been un- 
derestimated, in some cases substantially. On the other 
hand, in a dynamically active environment with poor 
shielding, an enhancement of VSGs relative BGs may be 
expected as enhanced erosion dust particle mantles turns 
BGs into VSGs. These arguments may also apply to the 
SMC. However, the 25/60 /im and 60/160 /im colors in 
NGC 1569 (« 0.13 and 1.2) are higher than in the SMC 
(« 0.05 and 0.4), more consistent with the hypothesis of 
a large population of very hot dust. Simultaneous fitting 
of the IRAC, MIPS, and IRAS data will be presented in 
a future paper. 

Resolved sub-mm data is cruci al to address th e most 
comp elling question raised by ILisenfeld et al.l ()2002l . 
l2005f ) : whe re, spatially, is the mi llimeter emission coming 
from? The lAguirre et al] (|2003l ) results, with a beam of 
several degrees, cannot address this question. At present, 
only one SMC molecular cloud has published, resolved 
mapping at millimeter wavelen gths. This is SMC Bl-1, 
mapped by SEST at 1.2 mm ()Rubio et all 120041 ). To- 
wards this cloud, we measure a background-subtr acted 
flux of Fbi_i,i6o ~ 1.5 Jy^ iRubio et all (|2004D find 
^Bi-1,1200 = 50 X 10-3 Jy. Therefore, i^i6o/^'i20o ~ 
30 for SMC Bl-1 compared to F160/-F1200 « 65 for 
the global S EP f see This CMC, described by 

iRubio et al.l (|2004f ) as cold and quiescent, shows a mil- 
limeter excess relative to the whole SMC, consistent with 
a cold dust temperature. However, the excess is to small 
to explain the the observed emission; extrapolating from 

SMC Bl-1 1% of the CO luminosity — to the 

whole SMC yields a millimeter flux of only a few Jy. 
lAguirre et all ()2003f ) measure 280 ± 80 Jy at A = 1.2mm. 

We present three explanations for the long wavelength 
emission: a wealth of cold dust, an abundance of hot 
small grains, and an emissivity law that differs from the 
Galactic one. We cannot rule out any of these, but prefer 
the latter, a modified emissivity, because it is simple and 

In a 1'.5 radius aperture abo ut CH950 = O'' 43™ 42^.4, (5ig5o = 
-73°35'10" l|Rubio et al.|[T993 ). with the background measured 
from a annulus between 5' and 6' to avoid contamination. 




Fig. 6. — The integrated SED of the SMC over the region shown 
in Figure|3]on hnear (top) and logarithmic (bottom) scales. We use 
DIRBE to fill in the (low intensity) regions be yond the coverage o f 
the MIPS maps. We scale the TopHat data lAguirre et al.ll2003l) 
by 0.9 to match the area we consider. The photometric consistency 
among the various bands is good. Nearby b ands agree within lo- 
in all cases. The solid line shows the best fit IDale fc H clou (200l) 
model, derived from the 100 and 160 ^tm fluxes. The three dotted 
lines show a single dust population with a temperature inferred 
from the lOO-to-160 /xm color and emissivity power laws with /3 = 
2, 1.5, and 1. 

because the available FIR/CO/mm data provide some 
evidence against the first and second explanations. 

3.4. Molecular Gas Map 

If dust and gas remain well mixed over the ~ 400 pc 
scales of our calibration regions, and the properties of 
grains do not vary too dramatically from the atomic to 
the molecular ISM, then the dust can be an optically 
thin tracer of molecular gas. Dust measured from FIR 
emission offers several key advantages over traditional 
molecular line tracers: unlike CO emission, dust emission 
is optically thin; unlike molecules, dust is not expected 
to be preferentially destroyed r elative to H2 in the oute r 
parts of molecular clouds (e.g. ' Malonev fc Blac3 ll988^: 
and the abundance of dust relative to gas may be directly 
measured elsewhere in the galaxy via comparison with 
the H I, while the use of molecular lines as tracers of gas 
must be calibrated within the GMCs t hemselves. In this 
section we u s e a m ethod adapted from 'Is rael (|1997[ ) and 
iDame et al.l ()200lD to construct a map of the H2 in the 
SMC. We summarize the method, present the map, and 
then discuss the uncertainties associated with the map. 

3.4.1. Method and Map 

To measure the H2 mass surface density, Sh2, along 
a line of sight we first adopt a dust-to-hydrogen ratio, 
D/H. We apply D/H to the dust mass surface density to 
estimate the total (atomic plus molecular) mass surface 
density along the line of sight, YiTotai- This is 



-'Total 



HI 



Sh2 — 



^Dust 

D/H 



(3) 



We remove the H I column, T^hi, from each line of sight 
to get a molecular gas column, Eh2 = ^Totai — ^hi- 
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TABLE 3 
Integrated SED of the SMC 
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Flux Density 
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HIRAS 


15500 ± 2200 
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IRIS 


14500 ± 2300 


100 


DIRBE 


15600 ± 2500 
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Literature" 


14000 - 20030 


140 


DIRBE 


16600 ± 4700 
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MIPS 


19000 ± 4300 


170 


Literature" 


15000 ± 2300 


170 


ISO 


18500 ± 4000 


240 


Literature" 


9600 - 12070 


240 


DIRBE 


10700 ± 3200 


476 


TopHat'' 


2800 ± 800 


652 


TopHat'' 


1400 ± 300 


750 


TopHat'' 


840 ± 200 


1224 


TopHat'' 


280 ± 80 


^ The ranf 


'e of values 


from the liter- 



ature compiled b y 
from ISchwerinp _ 
(|1988[). IStanimirovic et alj 



Wilke et a l 

IrIco ct al 

. (i2000l'i'. 

I Wilke et all (|2003t ). and lAguirre et al.l 



\2mc$ )}' Fluxes from TopHat scaled to 
he total fluxes from 
3200, 1620, 



amon area 



lAeuhre et all (1200311 
950, and 320 Jy. 



We m easure Sh/ from the map by IStanimirovic et alJ 
()1999f ). which includes zero spacing data. So 



^H2 



D/H 



(4) 



We measure D/H hy comparing the mass surface density 
of dust (FigureO to the mass surface density of H I over a 
series of reference regions, attempting to simultaneously 
satisfy two conditions: 1) that the reference regions is 
free of molecular hydrogen, and 2) that the reference 
measurement take place near the region where we want to 
determine Eh2) in order to minimize the effect of spatial 
variations in D/H. We adopt the following method: 

1. Wc identify regions likely to contain molecular gas 
from the CO map (Mizuno ct al. 2001, 2006). This 
is all lines of sight within 3'. 6 (65 pc) of Scr CO 
emission. We derive Sh2 over each of these regions. 

2. For each region, we create a reference region over 
which we measure D/H. This is all lines of sight 
more than 3'. 6 (65 pc) away from the CO but 
nearer than 10'. 8 (200 pc). The reference region 
is close enough that D/H should apply to the H2, 
but distant enough that contamination by H2 is 
unlikely. We check this by stacking CO spectra 
across all of our reference regions: we measure 
Ico ~ 0.04 K km s~^, so there is only very dif- 
fuse CO emission in the reference region (compare 
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Fig. 7. — This figure illustrates the methods used to construct 
the map of Eh2 shown in Figure |8] Thin black contours show 
CO emission (at 0.3, 0.6, 0.9, and 1.2 K km s~^), used to select 
regions over which we measure Sh2. Dotted contours show the 
boundary of the reference regions (merged together) over which 
we measure local values of D/H. We calculate Sh2 only within 
the inner dotted boundary; the reference region is free of H2 by 
construction. Crosses indicate selected peaks from both the CO 
and H2 maps that will be used below. The grayscale shows the 
derived Sh2. 

to Ico ~ 0.12 K km s"-"^ in the region outside of 
3cr CO but within 65 pc). 

3. Over each reference region, we measure the mean 
D/H and derive Eh2 using Equation U) 

Figure [7] illustrates this method. We show a map of the 
derived Sh2j the CO emission (black contours), and an 
outline of the set of all reference regions (dotted con- 
tours). Crosses mark peaks in the H2 and CO distribu- 
tions (used below). 

In order to match the resolution of the CO map, we de- 
rive Sh2 from a high resolution version of our dust map. 
We measure the lOO-to-160 /im ratio at the IRIS resolu- 
tion, but we derive the dust surface density (given the 
color) from the MIPS 160 ^m map at the CO resolution. 
We assume that the colors measured at 4' resolution ap- 
ply to the map at 2'. 6 resolution. We performed the same 
calculation using the HIRAS 100 /im map, which has a 
1'.7 resolution. We derive a comparable II2 mass, but the 
map has a number of artifacts with scales comparable to 
GMCs and therefore offers a check but no improvement. 

Figure [5] shows the molecular gas mass surface density, 
Sh2, derived using this method. The white contours in- 
dic ate the extent of 3cr signifi cant CO emission (0.3 K km 
s- MMizuno ern][2fml[200l assuming the average Xco 
we derive below (13 ± 1 x 10^^ cm~^ (K km s^^)^^). 
The II2 map contains a total of 3.2 x 10'' M0 of molec- 
ular hydrogen. We estimate the plausible range of II2 
masses given our systematic uncertainties as « 1.5 - 
4.0 X 10^ Mq. 

3.4.2. Uncertainties and Systematics 

The statistical uncertainty in Sh2 is small, only 
6 M0 pc~^. We estimate this from the measured uncer- 
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tainties in S_Dust and D/H and confirm it in two ways: 
1) by introducing noise at the 8% level (see Appendix 
[B|) into YiDust and deriving new Sh2 rnaps for many re- 
alizations, and 2) by computing I]h2 for several regions 
known to be free of CO emission. We also check the zero 
level in these empty regions, finding —5 to M© pc~^. 

Systematic uncertainties in the FIR maps — fore- 
ground, calibration, and saturation — will propagate to 
Sh2- To evaluate these, we offset the 100 /im and 160 
yum maps by ±1 MJy ster~^ one at a time and rederive 
Sh2- We perform similar tests to test the effect of our 
calibration and saturation; we scaling first the whole 160 
/xm map then just the bright regions by 1.25. Over all 
tests we find Mh2 in the range 1.4 x lO"^ to 3.9 x lO"^ Mq, 
suggesting a ^ 50% uncertainty. 

We avoid several systematic uncertainties by perform- 
ing the determination of the D /H locally for each region. 
Nonetheless, several astrophysical uncertainties may af- 
fect Eh2: the choice of dust map, changes in dust prop- 
erties between the atomic and molecular ISM, and the 
possibility of hidden cold dust. 

The uncertainty in Sh2 associated with our choice of a 
particular dust map is less than 30%. We estimate this 
using test maps adopting the various options described 
in Appendix |B] (i.e. single population, only the 60 and 
100 /im data, etc.). The total H2 inferred from resulting 
maps spans Mh2 = 2.1 - 3.2 x 10^ M©, with most val- 
ues between 2.4 and 3.2 x lO'^ M©. Note that applying 
this method directly to the MIPS 160 /Ltm intensity alone 
yields Mh2 = 3.8 x 10^ M©. 

Dust properties, especially the FIR opacity and D/H, 
may vary between the atomic and molecular ISM, pri- 
marily due to growth of icy mantles on the surfaces of 
dust grains in GMCs. These effects have the sense of 
decreasing I]h2 and a magnitude as high as a factor of 2 
- 3, though this must be considered an upper limit be- 
cause conditions in the SMC H2 will be less conducive 
to the growth of icy mantles than in Milky Way GMCs. 
See i)4.1l for more details. 

If the possible population of very cold dust discussed 
in §3.3.11 exists, our "Sdust map will be incomplete. If 
cold dust is concentrated in GMCs, we will underesti- 
mate T,£)ust and thus EH2 in these regions. In this case 
SH2 and S_Dust represent only lower limits. However, a 
large amount of cold dust is not required to explain the 
SED f ^3.3.ip and indeed we consider it unlikely. In ^4.2.21 
we find that quiescent gas with a high CO-to-FlR ratio 
represents « 30% of the CO emission in the SMC. If we 
allow that we miss all of H2 associated with these clouds 
in our map and assume that the average Xqo we mea- 
sure holds across the SMC, cold dust seems will represent 
a ~ 30% correction to the overall II2 mass. 

Table H] lists the uncertainties discussed here and gives 
our estimate of the maximum magnitude for each effect. 
UsTaeil (I199I gives a more thorough discussion of the un- 
certainties associated with this method. He concludes 
that most uncertainties have the effect that we underes- 
timate the true H2 mass. We noted two exceptions, in- 
creased FIR opacity and D/H in GMCs, above. Another 
is that we might mistakenly bookkeep warm ionized gas 
as H2. This requires very good resolution to investigate, 
so remains uncertain. H owever, studying the Eridanus 
Superbubble i Heiles et al.l (|T999i) found that IR not asso- 
ciated with H I is mainly associated with H2 , not ionized 
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+ 
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^ Per line of sight.'' The fraction of molecular gas cold 
enough that accretion onto grains is significant, which we 
argue to be small. 



gas. 

Several arguments suggest that our map of S]h2 is 
not dominated by these systematics. Most importantly, 
iDame et all (|2001f) found that H2 estimated from FIR 
and CO emission agree well in the Milky Way, so the 
method we apply to t he SMC work s in th e Galaxy. The 
estimate of Xqq that I Dame et al.l ()200lD derives using 
this method is 1.8 x 10^° cm~^ (K km s^^)~^, in ex- 
cellent agreement with estimates based on gamma rays 
fStrong & Mattox 1996). Further, Schn ee et all ()2006f ) 
studied the relation between FIR emission and dust col- 
umn for a simulated turbulent, externally heated cloud. 
They found that the 100/240 /J,m color (the best ana- 
log to our 100/160/im) does a good job of matching the 
extinction for Ay ^ 4, though it may underpredict the 
dust content above that value. Finally, Figure [8] shows 
good agreement between our H2 map and the CO map, 
particularly towards regions of active star formation. In 
fact, toward these regions the H2 is more extended than 
CO, suggesting that when dust is warm enough to emit 
in the FIR, the dust may trace H2 better than CO. 

Our estimate of the likely uncertainty in the I]h2 map 
is 50%. Systematic uncertainties in the data, dust mod- 
eling, and a modest population of cold dust are all im- 
portant but many of these effects have offsetting senses. 

4. ANALYSIS 

4.1. Dust and Atomic Gas 

Over the region of the MIPS 160 fj,m map, the SMC 
has an HI mass of Mhi = 1.8 x 10^ M©. The re- 
mainder of the 4 X 10^ M© of H I associated with 
the SMC (jStanimirovic et al.lll999f) lies in an extended 
component not covered by the MIPS map. With 
Mdust =3 X 10^ M©, this implies \og^^D/HI = -2.78, 
or l-to-600 over the whole SMC. If we include Mh2 ~ 
3 X 10^ then log^^ D/H = -2.86, or l-to-700. Includ- 
ing helium, the total dust-to-gas ratio is logj^Q DGR « 
-3.0. This is lower than the DGR in the Milky Way 
by about the same factor as the metallicity difference 
between the two systems with an uncertainty of a fac- 
tor of 2 (see Appendix |B|). This result is c onsis- 
tent with measurements by iBouchet et al.l ()1985l ). who 
fou nd a redden i ng per gas column 1/8 Galactic, and 
bv iWilke et all (|200l . who found a DGR of l-to-540. 
There are several predictions for how the D GR varies 
with Z. If we formulate DGR oc then llssa et al.l 
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Fig. 8. — The color image shows molecular gas mass surface density in the SMC inferred from the dust, from 10 to 200 Mq pc . The 
la uncertainty in Sh2 is o'eh2 ~ 6 Mq pc~^. The contours show the extent of significant (3(t) CO emission, which is 0.3 K km s~^, or 
» 62 Mq pc-2 at Xco = 13 ± 1 X 10^^ cm'^ (K km s-i)-l(see ^4.2. Ij l. Both maps are at the 2'. 6 resolution of the CO map. The dashed 
white contour shows the extent of the MIPS 160 fira map. 



(p9^ fin d a = 1.2, [Sc hmidt fc BoUeij (|l993l) fin d a = 
1.6, Lisenfeld fc Ferrara (.1998) find a = 1.9, and iDwekl 
(fl998[) find a = 1.3. These fits predict logio DGi? 
—3.0, —3.0, —3.3, —3.5. We cannot rule out any of these 
value s, but our derived PG R lead s us to prefer the shal- 
lower Issa et al. (1990) and lDweld ()1998[ ) fits. 

Figures [5] and 1101 compare dust and atomic gas in our 
maps. Figure [10] shows that dust surface density and 
atomic gas surface density are correlated across the SMC, 
with a rank correlation coefficient ^ 0.76 across all 
data. All three figures show that the D/HI is a function 
of location in the SMC (previously n oted by several au- 
thors and shown nicely in profile by IStanimirovic et al.l 
[2000l) . Diffuse gas in the Wing — distinguished by a 
simple intensity cut, /igo > 10 MJy ster~^ (which is 
arbitrary but does a good job) — shows a low mean 
D/HIwing = -2.94 ± 0.02, while for the Bar 
\ogiQ D / Hlsar = —2.71 ± 0.01, a factor of two higher. 



Lines of sight towards regions with CO emission show 
higher D/HI than nearby lines of sight, the mean is 
logio D/HIco = -2.54 ±0.02. It is therefore possible to 
identify most of the SMC star-forming regions from the 
map of D/HI, Figure [1 

Some of the variations in D/HI r nay be understood 
in terms of the life cycle of dust. iDwelJ (|T998;1 out- 
lines a picture in which dust is depleted in super- 
nova sh ocks and accrete s new material within molecular 
clouds. iBot et all ((2004) argued that supernovae shocks 
will destroy dust more easily in the diffuse gas found in 
the Wing. Semb a ch fc Savage (1996) observed a similar 
effect in the Milky Way, that depletion of some elements 
decreases with increasing distance from the plane of the 
Galaxy. This effect, combined with the distance from the 
sources of new dust in the Bar, may explain the enhanced 
depletion in the Wing. 

The values of the D/HI near CO emission must be 
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due, at least in part, to the presence of molecular gas; 
the presence of CO emission necessitates some H2 along 
the line of sight. If this is the dominant effect, then H2 
accounts for ^ 30% of the gas surface density in these 
regions, on average. The processing of dust inside GMCs 
may also contribute to the high logj^g D /HIco ■ accretion 
of heavy elements onto grains may result in an enhanced 
d ust abu ndance and higher FIR opacities. 

iDwekl ([l998l estimates that over the lifetime of a Milky 
Way GMC, the dust mass in the cloud may increase by 
a factor of 2 - 3 due to accretion of heavy elements (see 
also references therein); so for Milky Way GMC in the 
middle of its lifetime, we would expect an enhancement 
in the DGR of ~ 1.5 - 2. However, an observational 
measure of the average DGR enhancement between the 
diffuse and molecular ISM is har d to come by. In the 
Milky Wav. lRachford et al.l ()2002D found gas and dust to 
be well mixed along lines with Ay ^ 1. There is some 
evidence from absorption studies that the DGR varies in 
the SMC with higher values toward lines of sight with 
more II2 ( Sofia ct al. 2006), but the sample is too small 
to evaluate the significance of this finding. 

Some of the enhanced D/HIco may also result from 
dust in GMCs having larger FIR emission per unit 
mass at a given temperature (i.e., higher kfir). In the 
Milky Way, dust associated with dense, colder gas is ob- 
served to have a higher emissivity than diffu se, warm gas 
(jCambresv et al.l 120011: IStepnik et all 120031) . Modeling 
the effect of icy mant l es on grains in protostellar cores, 
lOssenkop f fc Henningl (|1994l ) found an enhancement of a 
factor of 2-3 in the FIR opacity. 

Both of these effects depend on conditions found at 
high Ay which may be more common in Milky Way 
GMCs than in the SMC. We shall see below that Ay 1- 
2 is typical for an SMC GMC, ^ 3 times less on average 
than for Milky Way GMCs. Therefore these conditions 
will not be as prevalent in SMC II2 as in the Milky Way. 
If DGR enhancements and increased kfir are important 
effects, then the true Ay towards SMC GMCs is even 
lower. Therefore, we argue that the factors of 2-3 dis- 
cussed here for both ef fects must be viewed as upper lim- 
its. In particular, the lOssenkopf fc Hennind ()1994l ) cal- 
culations are for cores with densities > 10^ cm~'^, which 
will represent only a tiny fraction of the molecular gas in 
the SMC. While all three effects contribute to the high 
value of D/HIco, we consider it likely that presence of 
II2 is primarily responsible for the high D/HIco- 

4.2. CO and H2 in the SMC 
4.2.1. The CO-to-H2 Conversion Factor 

Our II2 map contains 3.2 x 10^ Mfri of molecu l ar hy- 
droge n. The CO luminosity in the iMizuno et"an (|2001l . 
|2006[) NANTEN map over the same area is 1.5 x 10^ K km 
s~^ pc~^ (summing over the entire map without mask- 
ing). This implies a mean C0-to-H2 conversion factor, 
Xco, of 13 ± 1 X 10^1 cm-2 (K km s"i)"\ which is 65 
times the Galactic value. This conversion factor varies 
systematically within the SMC. Regions with high CO 
intensity [Ico > 0.3 K km s'^) have Xco = 9± 1 x lO^^ 
cm~^ (K km s~^)~^, while regions with diffuse CO emis- 
sion [Ico < 0.3 K km s~^) have twice this value, 
Xco = 18 ± 2 X 10^1 cm-2 (K km s"^)-^. 

If CO and II2 are not coincident (see i )4.2.2|) . then 
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Fig. 9. — The ratio of dust mass surface density to atomic gas 
mass surface density at 4' resolution. Star-forming regions show 
very high dust-to-gas ratios, probably because they contain molec- 
ular gas, which is not accounted for in this map. The star-forming 
Bar of the SMC shows higher D/HI than the diffuse Wing. 
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Fig. 10. — Dust surface density as a function of H I surface den- 
sity. Lines of sight towards CO emission, the SMC Bar, and the 
low intensity Wing are marked with black, light gray, and dark 
gray dots respectively. Lines of constant dust-to-gas ratio are la- 
beled. The correlation between dust surface density and H I is good 
across the SMC. The Bar shows a higher D/HI than the Wing, 
and regions near CO emission show the highest dust-to-gas ratios, 
probably because they harbor significant amounts of molecular gas. 

we must be careful in defining exactly what is meant 
by the CO-to-Il2 conversion factor. Indeed, much of 
the controversy in the literature is caused by confus- 
ing large-spatial scale conversion factors, which may in- 
clude II2 without associated CO, with local conversion 
factors, which are limited to regions with substantial 
CO emission. Therefore, we calculate a refined value 
of Xco over the common volume shared by CO and II2 
for the molecular peaks marked in Figure [71 We com- 
pare CO emission to H2 in a 1' diameter aperture, es- 
sentially measuring the beam pointed at the center of 
the cloud. We then adjust the H2 down by a factor of 
1.3 to account for H2 along the line of sight that lies 
outside of the CO emitting region (we find in i )4.2.2l 
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that the H2 is more extended than the CO by this 
factor). Table [6] shows the results. We find a mean 
Xco = 6 ± 1 X 10^1 cm-2 (K km s-^)-\ We derive a 
similar value, Xco = 7 ± 1 x 10^^ cm^^ (K km s~^)~^, 
by considering lines of sight with bright CO {Ico > 0.6 K 
km s""'^) and scaling by the radius ratio (1.3). Because 
it represents our best estimate of Xco over the volume 
from which the CO emission originates, this is the value 
most directly comparable to virial mass estimates. 

Table [5] summarizes Xco from this paper and the lit- 
erature. Most of these values are from virial mass mea- 
surements, which obtain values of Xco that are sys- 
tematically higher than the Galactic value but depen- 
dent on the scale of the structure studied (jRubio et al.l 
ll993tlBolatto et aLll2003f ). Larger structu res yield higher 
values of Xco- Indeed, iBolatto et al.l ^2003), study- 
ing the N 83/N 84 region, found nearly Galactic values 
of Xco when considering the smallest resolved struc- 
tures in their data set (GMCs with radii of 10 to 20 
pc), but at the scale of the entire complex found a con- 
version factor ~ 10 times the Galactic value. At a 
common scale of ~ 40-100 pc, the virial mass method 
yields Xco = 0.9-3.5 x lO^i cm-^ (K km s'^)-^. At 
the same scales, our FIR based II2 map yields Xco ~ 
6 ± 1 x 10^^ cm~^ (K km s~^)~-^, 2-6 times higher. 

The discrepancy between Mh2 measured from the FIR 
and the virial mass may be explained if magnetic fields 
may provide most of the support for SMC GMCs. Mag- 
netic support will not affect the turbulent line width, 
so in this case the virial ma ss will und erestimate the 
true mass of th e CMC (e.g. iMcKee fc Zwcibcl 1992). 
iBot et all (l2006h compare H2 masses estimated from mil- 
limeter emission to CO-based virial masses for a set 
of SMC and Milky Way GMCs. They find that in 
Milky Way clouds the virial mass usually exceeds the 
millimeter-based mass, but that SMC clouds have higher 
millimeter masses than virial masses. They attribute this 
difference to enhanced support by magnetic fields in SMC 
clouds. If magnetic support does not play an important 
role, SMC clouds may be short-lived as this measurement 
implies that the kinetic energy is much smaller than the 
g ravitat io nal po tential energy. 

[Israeli ()1997[ ) applied a method similar to ours 
to estimate Mh2 towards several regions in the 
SMC with a limiting resolution of 15' (^ 270 pc). 
Given this resolution, his value of Xco= 12 ± 2 x 
10^^ cm~^ (K km s~^)~^is best compared to our global 
value, 13± 1 x 10^^ cm"^ (K km s~^)~\ and the two 
are identical within the uncertainties. Application of 
his method to the newer data used in this paper yields 
Xco = 9 - 13x 10^1 cm-2 (K km s-^)-^ for lines of sight 
with CO emisqsion, depending on the method used to de- 
fine reference regions. The similarity betwe en our results 
using the 100 and 160 /xm bands and the llsraell ()1997l ) 
results using the IRAS 60 and 100 /im bands shows that 
a color-corrected 100 fim map is enough to capture the 
gross behavior of H2 in a galaxy. This is also seen by 
the success of lDame et al.l (j2001.) in matching Galactic 
FIR emission to the CO distribution. Longer wavelength 
(especially sub-mm) data may improve the accuracy of 
column density es timates by ~ 10% a nd reduce the scat- 
ter substantially (|Schnee et al.l 12006') . but is not neces- 
sary to obtain the basic result. The improvements in the 
T,H2 map presented here are chiefly: resolution capable 



TABLE 5 
CO-TO-H2 Conversion Factor 



Source 






cm^'^ (K km s"-"-)^^ 


Whole SMC'' 


13 ± 1 X 10^1 


LOS Ico > 0.3 K km s-^^^ 


9 ± 1 X 10^1 


LOS Ico < 0.3 K km s-^'' 


18 ± 2 X 10^1 


Molecular Peaks^'*' 


6 ± 1 X 10^1 

1 c in2i 




Individual GMCs 
... N83/N84 Complex 
Israel figQ?) "-" 
Rubio et al. ('20041° 



0.9-1.5 X 10^1 
9.0 X 10^0 (_R/10 pc)°-'' 

1.8-7.8 X 10^0 
2.2 X 10^1 



1.2 X 10^^ 
1.3 - 3.6 X 10^2 



Using FIR to trace H2.'' Marked in Figure [7] Ad- 
justed by 1.3 to account for extended H2 beyond 
CO along the line of sight. Most appropriate for 
comparison with the virial mass.'^ Using virial mass 
method.^ Using IRAS 60 and 100 ^im bands." Sub- 
millimeter continuum observations. 

of distinguishing Galactic GMCs (~ 50 pc), allowing us 
to study the relative structure of CO and H2; and the 
use of local reference regions to calibrate out vari ations 
in D/ H across the SMC — the low resolution of Israeli 
()1997t ) made a single D/H appropriate for that work. 

4.2.2. Comparison Between H2 and CO 

Figure ID shows that CO predicts the presence of II2 
almost perfectly. However, the reverse is not true; we 
derive substantial II2 surface densities beyond the 3a CO 
emission (white contours). Thus, H2 appears more ex- 
tended than CO emission in SMC star-forming regions. 
However, an alternative explanation is hard to rule out: 
because Xco is very high, the surface density sensitivity 
of the CO map may in fact be quite poor. In this case, 
there may be low intensity CO below the noise level of 
the map and the differences in the area are an artifact of 
sensitivity. In this section we show that these differences 
are more than a sensitivity effect; there are substantial 
variations in Xco across the SMC. 

Figure [T2| shows the distribution of H2 and CO about 
several molecular gas peaks. These peaks are indicated 
in Figure |7| and summarized in Table |6l We normal- 
ize each profile by its peak value, so Figure [T^j contains 
no information about the relative amplitude of the two 
distributions, i.e. Xco- Instead, Figure [T2l shows the 
relative structure in the two maps. For four of the six 
peaks, H2 is more extended than CO; around one peak 
CO is more extended than H2; and around one peak, the 
N 83/N 84 region in the east, CO and H2 have similar 
structure. 

We stack the normalized profiles and measuring the 
half-width at half maximum (HWHM) for the average 
H2 and CO profiles. The stacked H2 has HWHM- 41 pc 
while the stacked CO has HWHM'^ 30 pc. Integrating 
the average H2 profile and the average CO profile we find 
/ SH2/ / ^co — 1-3- Thus, along an average Hne of sight 
towards a molecular peak, — 30% of the H2 lies outside 
of the region with CO emission — a number we used to 
compute the "common volume" Xco estimate above. 

Figure [11] shows another view of the relative distri- 
bution of H2 and CO. We plot regions with Sh2 > 



16 



Leroy et al. 



60 M0 pc^^ in gray and indicate the extent of signifi- 
cant CO emission (Scr = 0.3 K km s^^) with black con- 
tours. These two thresholds arc matched for Xco = 
13 ± 1 X 10^^ cm~^ (K km s^^)^^, our integrated con- 
version factor. Note that this is ^ 2 times higher than 
our estimate of Xqo over the shared volume. Therefore, 
the gray area in Figure [11] shows a very conservative es- 
timate of where would expect to find CO emission based 
on our FIR-derived Eh2- Despite this conservative cut, 
toward regions of active star formation — at the north 
and south of the Bar and in the eastern N 83 region — 
H2 extends beyond the CO. In these regions, the dust is 
warm and Figures [Tl] and [8] give us confidence that our 
FIR-based H2 map traces SH2 better than CO. 

Thus we present three measurements that suggest the 
H2 to be more extended than the CO: 1) the profiles 
shown in Figure [T^l 2) the large extent of CO relative 
to H2 at a common, conservative threshold in Figure [TTl 
and 3) most compellingly, the measurement of different 
values of Xco averaging over low and high Ico regions 
— where /co > 0.3 K km s'^ we find Xco = OilxlO^i 
cm~^ (K km s~^)~^ and where Ico < 0-3 K km s~^ we 
find Xco = 18 ± 2 X 10^1 cm'^ (K km s'^)-^. None of 
these measurements depend critically on the amplitude 
of the H2 map, i.e. Mh2, or the sensitivity of the CO 
map, but only on the relative structure of the H2 and 
CO maps. We interpret this increase in H2 relative to 
CO towards the outer parts of SMC clouds as evidence 
of the selective photodis sociation of CO in these regions 
iMalonev fc BlaS (flQSl . 

Away from the star-forming regions, there is more CO 
relative to FIR-derived II2 . Several small CO clouds — in 
the Wing, the central Bar, and outside the northern and 
southern star-forming regions — may be good candidates 
to harbor populations of cold dust. In these regions, we 
still find significant II2 emission, — 90% of lines of sight 
with Ico > 0.3 K km s~^ also have T,h2 > 20 Mq pc"^ 
(3a) and the remainder all have Eh2 at slightly lower 
significant. However, towards these regions we may un- 
derestimate the total Sh2 because there is a population 
of cold dust that is incompletely traced by our FIR maps. 

How much molecular gas is represented by these qui- 
escent clouds? The region with I]h2 > 60 M0 pc~^ con- 
tains both 70% of the significant CO emission and 70% 
of the significant H2 surface densities. While the FIR- 
based I]h2 may fail in these regions, they will be subject 
to less intense radiation fields and should therefore suffer 
less selective photodissociation of CO. Therefore, based 
on the fraction of the total CO luminosity coming from 
such clouds, we estimate the fraction of the SMC molec- 
ular phase in cold, quiescent clouds to be ^ 30%. 

4.2.3. Surface Density and Extinction 

If photoionization by far-UV photons regulates the 
structure of molecular clouds, then GMCs in low metal- 
licity systems may be expected to have higher surface 
densities than Milky Way GMCs in order to achieve 
the same mean extinction as their Galactic coun terparts 
(average Ay sa 4 - 8 over a Milky Way CM C, iMcKeel 
119891: ISolomon et 311119871 ). iPak et al.l (il998D estimated 




Fig. 11. — The extent of CO and H2 from the FIR at a com- 
mon significance assuming our global CO-to-H2 conversion fac- 
tor, Xco = 13± 1 X 10^^ cm-2 (K km s"i)-i. Both maps are 
clipped at the 3a for the CO, or Sh2 ~ 60 Mq pc~^. In regions 
of active star formation, the H2 shows a greater extent than the 
CO, but away from these regions, CO emission may indicate the 
presence of cold(er) dust. 

We find surface densities similar to those of spiral 
galaxy GMCs rather than th e val ues of > 1000 pc~^ 
expected from|McKed (|1989( ) and lPak et all (|1998f i. The 
molecular gas peaks indicated in Figure [7] have surface 
densities I;h2 = 90 - 310 with a mean Eh2 = 180±30 Mq 
pc~^. This similar to the mean for Mi l ky W ay GMCs, 
Sh2 ~ 170 Mq pc-2 (So lomon et al.lfl98l . If CO- 
dark H2 ma kes up 10 - 50% of the molecular gas in the 
Milky Way (jGrenier et al.ll2005f ). SMC clouds may even 
have 'Sh2 slightly lower than Galactic clouds. Along the 
same lines of sight, the H I surface density varies from 
= 35 - 105 Mq pc~2 with a mean of E/// « 70 ± 10 
Mq pc~^. Over the entire region with 3a CO emission — 
a reasonable estimate of the star-forming area — we find 
a mean I]h2 = 90±10 Mq pc~^ and a mean Y.hi = 65±2 

Mq pC-2. 

Nearly Galactic surface densities imply that the mean 
extinction through SMC GMCs i s much lower than in 
Galactic GMCs. iBouchet et all (|1985[ ) found Ry = 
Av/E{B - V) = 2.7 ± 0.2 for the SMC. Using this 
value and Nh/E(B - V) = 5.8 x lO^^ cm^^ niag-^ 
(|Bohlin et al.lll978f) . scaled down by the D/H relative 
to Galactic (taken here to be « 0.01), we calculate Ay: 



A, 



NjHI) + 2 N{H2) D/HsMC 
2.1 X 1021 cm-2 D/Hmw 



(5) 



N{H) > 1023 cm-2 (5^^ > goo Mq pc-2) towards SMC 
star-forming regions. Using our dust and gas maps, we 
calculate the mean surface density and extinction to- 
wards the molecular peaks in the SMC. 



We measure D/Hsmc/D/Hmw to be « 7 (322]). Us- 
ing Equation [5] we find Ay ~ 1.8 ± 0.6 mag towards 
the molecular peaks and Ay « 1.1 ± 0.5 mag over the 
w hole region wi t h CO emission. Apply the method 
of iSchlegel et all (|1998D to our dust map yields consis- 
tent extinctions, Ay r:! 1.6 ± 0.3 toward the peaks and 
Ay w 1.0 ± 0.4 over all CO emission. 

Our H2 map has 46 pc resolutio n, close to the averag e 
diameter of a Milky Way CMC ()Solomon et al.lll987^ . 
If the star-forming parts of SMC GMCs have the same 
diameters as Milky Way clouds, then our data are not 
consiste nt with the high surface densities predicted by 
iMcKed ([1989, ) for low metallicity GMCs. However, the 
region of CO emission, rather than the more extended H2 
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TABLE 6 

Properties Of Several Molecular Peaks 



a 

(J2000) 


5 

(J2000) 


^HI 

(Mq pc-2) 


^H2 

(M0 pc-2) 


lO^^cm ^ 
K km 


HWHM(H2) 
(pc) 


HWHM(CO) 
(pc) 


logio D/H 
(Reference) 


01hl3m45.8s 


-73dl6m30s 


58 


261 


7.1 


40 


36 


-2.93 


01h02m49.3s 


-72d02m51s 


49 


133 


10 


44 


33 


-2.70 


00h58m54.8s 


-72d09m24s 


36 


81 


5.4 


35 


22 


-2.64 


00h48ml0.4s 


-73d06m03s 


106 


230 


5.6 


53 


22 


-2.82 


00h47m46.4s 


-73dl5m34s 


83 


162 


4.0 


69 


26 


-2.82 


00h45m26.2s 


-73d22m45s 


70 


199 


6.5 


40 


54 


-2.82 


^ Applying a 


correction of 1/1.3 to account for H2 


beyond the CO along the line of sight (see 


text). 
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Fig. 12. — Normalized CO intensity (black) and H2 surface density (gray) in concentric rings around selected maxima. All profiles are 
normalized to the central value and then averaged in concentric rings. The dashed lines show the profile of a point source. The H2 profile 
derived from the FIR and H I is more extended than the CO emission in 4 of 6 cases; in the other two cases the profiles are comparable. 
The position of each maximum is indicated above the panel and the maxima are noted with green crosses in Figure |8] 
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region, may be the relevant comparison for the iMcKed 
(|1989D prediction. We only expect star formation to oc- 
cur within the CO-emit t ing r egion, and the Ay = 4 - 
8 calculated by iMcKed (|l989f ) is for GMCs defined by 
CO emission. In this case, resolution effects ("beam dilu- 
tion" ) may cause us to underestimate the surface density 
associated with the CO peaks. If the radius of a typical 
SMC CMC is ~ 15 pc, for example, then the associated 
surface densi ty of may be ^ 1600 Mq pc~^, consis- 
tent with the'McKee fc Zweibell jl992) prediction. Many 
structures in the Mizuno et al.l ([2006.) CO map are un- 
resolved or only marginally resolved, and other studies 
have indicated that CO clouds in both LMC and SMC 
may be smaller than in the Milky Way, with CO confined 
to small re gions and little o r no diff use component ( e.g. 
iLeaueux et al...l994l : lRubio| [l999: Mizuno et al.|[2001l) . 

Thus, we find evidence for differences between SMC 
and Milky Way GMCs, but cannot distinguish the sense 
of them because of limited resolution. SMC GMCs either 
have lower average extinction through the clouds or are 
substantially (^ 3 times) smaller on average than Milky 
Way GMCs. A second high resolution band is needed to 
test which of these conclusions holds by measuring the 
surface densities of SMC GMCs at a resolution of « 10 pc 
(that of the MIPS 160 ^m map). 

4.2.4. Pressure and Molecular to Atomic Gas Ratio 

In spiral galaxies there is a good correlation between 
the ratio of molecular gas mass to atomic gas mass and 
the h ydrostatic midplane pressure (|Blitz fc Rosolowskvl 
l2004f ). Dust is integral to H2 formation and may impede 
H2 destruction (though self-shielding of II2 may also have 
this effect). Therefore we expect that, given the low dust- 
to-gas ratio in the SMC, the ratio of molecular gas to 
atomic gas for a given pressure (gas density) will be lower 
i n the SMC than in sp iral galaxies. 

iBlitz fc Rosolowskvl (j2004) derive that the hydrostatic 
midplane pressure for a thin disk of gas in a larger disk 
of stars is 

= 0.84(GS,)"-5Sg^^ (6) 

where Eg is the total surface density of the g as , Xj ^ IS 
the surface density of stars, Vg is the velocity disper- 
sion of the gas, and ft,* is the scale height of the stellar 
disk. The applicability of this formula in the SMC is 
questionable: the three dimensional distribution of the 
stella r potential is l ikely not well described as a disk. In- 
deed, ICrowl et al.l ()200l f) find a la dispersion of 6 - 12 
kpc along the line of sight through the SMC. They sug- 
gest that the galaxy may be a triaxial spheroid with a 
1:2:4 ratio of sizes along the right ascension, declination, 
and line-of-sight directions. Also, although there is usu- 
ally only a single molecular complex along a given line 
of sight, there are often several H I components. Even if 
Equation [6] applies, the inclination, stellar scale height, 
and the mass-to-light ratio of the stellar and gas disks 
are not known to better than a factor of two. 

Emphasizing the approximate nature of the calcula- 
tion, we derive the pressure and the ratio of molecular 
to atomic gas at ^ 0°.7 resolution 750 pc). The 
mean velocity dis persion of H I along the li ne of sight is 
Vg « 20 km s"i (jStanimirovic et al.|[l999l ). We adopt 



ft* = 6 kpc (the low end of thelCrowl et al.|[200lL range) 
and no inclination. We use the DIRBE if-band map 
and a mass-to-ligh t ratio M/Lk — 0.5 Mq/Lk.q (e.g. 
ISimon et al.l l2005l ) to calculate E*, and the FIR-based 
II2 map combined with the H I map to calculation the 
gas surface density and the ratio of molecular to atomic 
gas. 

Figure[T3]shows fmoi as function o f Ph for the SMC and 
a samp le of spiral galaxies studied bv lBlitz fc Rosolowskvl 
()2006f ) with similar spatial resolution. Most lines of 
sight for the SMC have Ph/kB = 10^ - 10^ K cm'^, 
with almost half of the data in the bin centered on 
PhlkB = 2.4 X 10"! K cm -3. In this bin, the best fit 
iBlitz fc Rosolowskvl (|2006D relation predicts fmoi ~ 0.6 
for large spirals. For reference, this Ph and fmoi are sim- 
ilar to those found at the Solar Circle in the Milky Way, 
a region that is dominated by H I by about 2-to-l. In 
the SMC we measure fmoi = 0.2 in the bin centered on 
Ph/ks = 2.4 X 10^ K cm~'^, a factor of 3 lower than the 
prediction/Milky Way value. Over the whole SMC, we 
find fmoi a factor of 2 - 3 lower than predicted. 

Thus, Figure [T3| suggests that the low ratio of H2 to 
H I in the SMC is due to a combination of low gas den- 
sities and environmental effects. Because the median Ph 
in the SMC is lower than in the inner disks of spiral 
galaxies, fmoi is low. However, fmoi is lower than found 
in spiral galaxies by an additional factor of 2 - 3. This 
may be the result of the low DGR. We suggest the fol- 
lowing explanation: Because the formation rate of II2 
is limited by the rate at which H ato ms stick to the 
surface of dust grains (HoUcnbach. Werner, fc Salpeteij 
|1971| ). II2 will form more slowly in the dust-poor SMC. 
This lower formation rate means that given the same 
conditions — density, shielding from dissociating radia- 
tion — the equilibrium II2 /H I ratio should be lower for 
a low DGR galaxy. Alternatively, SMC clouds may be 
out of equilibrium: for a low DGR, an H I cloud will take 
longer to convert to H2, even under conditions where an 
overwhelmingly molecular cloud is the equilibrium state 
(jBell et al.ll2006t ). If the timescale to reach equilibrium 
is long compared to the cloud lifetime, then these clouds 
may never reach chemical equilibrium. In this case, we 
may expect to find a lower molecular gas to atomic gas 
ratio in low mctallicity galaxies like the SMC. 

5. CONCLUSIONS 

We present new FIR maps of the SMC at 24, 70, and 
160 /xm obtained by Spitzer as part of the S'^MC survey 
(|Bolatto et al.ll2006D . These maps cover the Bar and the 
Wing regions with 6", 18", and 40" resolution at 24, 70, 
and 160 ^m respectively (1.8, 5.4, and 12 pc). Combining 
these maps with data from the literature, we draw the 
following conclusions: 

1. We find a total Mdust =3 x 10^ Mq, which is sys- 
tematically uncertain by a factor of 2 — 3 because 
of degeneracies in how to model the dust. Strong 
upper and lower limits to the dust mass in the SMC 
are ~ 2 x 10'' and 10^ M©, respectively. 

2. If the long wavelength part of the SMC SED comes 
from a population of very cold dust, then Mdust is 
a lower limit. However, based on the lone GMC 
with published, resolved millimeter observations 
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Fig. 13. — The molecular to atomic gas ratio, /moli as a function 
of hydrostatic midplane pressure, conducted at 0° .7 resolution. We 
use the formula of Blitz & Rosolowsky (2004) with o-„ = 20 km s~^ 
and h* = 6 kpc and also plot data for spiral galaxies measured 
at comparable spatial resolution (Blitz & Rosolowskv 2006). The 
histogram above the plot shows the fraction of lines of sight in each 
pressure bin for the SMC. The low H2 to H I ratio in the SMC is 
a result of the relatively low values of and an additional factor 
of ~ 3 deficit that may be the result of a low dust-to-gas ratio. 

it seems unlikely that GMCs harbor the requisite 
amount of cold dust to account for the long wave- 
length emission. There are plausible alternative ex- 
planations: an abundance of hot VSGS or that the 
emissivity of SMC grains at long wavelengths dif- 
fers from the Galactic value. The latter explana- 
tion is simple and consistent with the data. 

3. The dust-to-hydrogen ratio, D/H, over the whole 
region studied is log^Q D/H — —2.86, or l-to-700, 
including H2, implying a total log^o DGR w —3.0. 
This is about Zsmc/Zmw times the Milky Way 
PG R and is consistent with iBouchet et al.l ()1985D 
and lWilke et"aI1 (I200I . 

4. The dust-to-H I ratio, D/HI, is a function of po- 
sition in the SMC. In the Wing and the Bar, we 
find mean log^Q D / Hiding = —2.94 ± 0.02 and 
logio D/iJ/sar = -2.71±0.01. D/ff/ near regions 
with CO emission is higher still: \ogiQ D/HIco = 
—2.54 ± 0.02. Some of this high value may be due 
to enhancements in the opacity and abundance of 
dust in molecular clouds, e.g. due to growth of 
icy mantles. However, we argue that most of the 
D/HI enhancement near CO emission is best ex- 
plained by the presence of molecular hydrogen. 

5. Following a i netho d outhned by llsraell ()1997f ) and 
iDame et al.l ()2001f ). we derived a map of I]h2 in 
the SMC. To derive this map, we calibrate D/H 
locally using a set of reference regions show in 
Figure [7] and apply Equations [H] and [H We find 
AfH2 = 3.2 X 10^ Mq of H2 in the SMC, equal to 
~ 10% of the H I mass. This is lower than the ra- 
tio of molecular to atomic gas found in the inner 
parts of spiral galaxies, but similar to that found 



i n dwarf galaxies and the outer parts of spirals 
(lYoung fc Scovilld flMt Hsraell [l997l : iLerov et all 
12005!, and see Figure [HI). 

6. The derived C0-to-H2 conversion factor over the 
entire SMC is 13 ± 1 x 10^ | cm- ^ (K k m s'^)-^, 
similar to that obtained by [Israeli (|1997D . Toward 
molecular peaks, correcting for II2 beyond the CO 
along the line of sight, we find a conversion factor of 
6 ± 1 X 10^^ cm~^ (K km s~^)~^ This represents 
our best value of Xqq for comparison with virial 
mass measures. It is still 2-4 times larger than 
virial Xco measurements for SMC GMCs. 

7. SMC GMCs appear to be gravitationally bound. 
We measure Mh2 to be larger than M^^^, surface 
densities similar to those of Milky Way GMCs, and 
a mean midplane pressure to similar to that at the 
Solar Circle (albeit with a large systematic uncer- 
tainty) . 

8. We find evidence for the selective photodissociation 
of CO. About the molecular peaks, II2 is more ex- 
tended than CO by about 30%. Averaging over low 
CO {Ico < 0.3 K km s'l) and high CO {Ico > 0.3 
K km s~^) regions, we find the low CO regions have 
twice as much H2 per CO as the high CO regions. 
We average over a large region to make these mea- 
surements, so this is not a sensitivity effect. Re- 
gions of high H2-to-CO occur near intense radia- 
tion fields surrounding bright CO clouds, consistent 
with the scenario in which the outer parts of SMC 
clouds see their CO destroyed by UV photons while 
II2 survives due to self-shielding (jMalonev fc Blacl3 
IT988) . 

9. Several clouds away from massive star formation 
show more CO relative to FIR than average, mak- 
ing them good candidates to harbor populations 
of cold dust. On e of these clouds, SMC Bl-1 
(jRubio et al.ll200^ , shows a greater submillimeter- 
to-FIR ratio than the SMC as a whole. These 
clouds represent ~ 30% of the total CO luminosity. 

10. GMCs in the SMC are either translucent or very 
small. SMC molecular peaks have H2 surface den- 
sities similar to Milky Way GMCs, I;h2 ~ 100 - 
200 M0 pc~^ and extinctions of Ay « 1 - 2 mag. 
This is lower than lower tha n the Ay = 4 - 8 mag 
found in Milky Way clouds ()McKeelll98"9l ). Either 
photoionization by far-UV photons does not regu- 
late the structure of SMC GMCs and they indeed 
have lower average extinctions or the star-forming 
parts of SMC GMCs are a factor of ~ 3 smaller 
than our 46 pc beam. 

11. The low Ha-to-HI ratio in the SMC results from 
a combination of low gas densities and metallicity 
effects. The SMC has a lower ratio of molecular to 
atomic gas by a factor 2-3 th a.n expected based on 
its h ydrostatic pressure (and iBlitz fc Rosolowskvl 
I2OO6 ), but also a median pressure a factor of 3 
lower than large galaxies. The median Ph is com- 
parable to that at the Solar Circle, though very 
uncertain. 
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APPENDIX 
CALIBRATING THE MIPS 160 MAP 

Figure[6]and Table[3]show that there is good photometric consistency among our data. Fluxes at the same wavelength 
agree within the uncertainties, and the data is consistent with DIRBE, which represents the best measurements of 
absolute sky brightness at FIR wavelengths. The exception to the general good agreement is that the 160 /im contains 
more flux than we expect based on comparison with the DIRBE data at 140 and 240 and the ISO 170 /im 
map. In this section we describe how we use the DIRBE data to calibrate the MIPS 160 /im map. The calibrated 160 
Hm map is used for all of the results in this paper, including Figure [5] and Table [3] Our approach to calibrating the 
160 fiia map is: 

1. We adopt DIRBE as our absolute calibration. 

2. From the DIRBE bands at 140 fim and 240 /im we estimate the intensity at 160 /im by interpolation. We adopt 
(3 = 2 and calculate a temperature from the ratio of the 140 /im to 240 /im bands. The method of interpolation 
makes little (< 10) difference because the 140 /im DIRBE band is close to the 160 /im MIPS band. 

3. The foreground-subtracted DIRBE data (see Sj2]) are our estimate the 160/im intensity from the SMC along each 
line of sight. We add our best estimate of the MIPS 160 /im foreground (from zodiacal light, Galactic cirrus, 
and the CIB) to this value. This is our best estimate of the intensity that MIPS should measure. 

4. We make the decision to model the difference between MIPS and DIRBE using one parameter: the best fit gain. 
This is the ratio of the intensities measured by MIPS to those that it should measured as estimated from DIRBE, 
Idirbe I Imips 

5. We derive the best value for the gain by taking the average ratio of predicted to measured intensity over the 
MIPS map. We perform the comparison at the DIRBE resolution, 0°.7. We use DIRBE to fill in empty parts of 
the MIPS map, but we only considering lines of sight for which at least 75% of the data within the beam come 
from MIPS. This average ratio is 0.8, so that lDiRBE,i60tim — 0.8 Imips, i60fj.m', the median and integrated ratios 
within a few percent of this value. 

6. We apply this factor to the MIPS data before foreground subtraction and proceed as described above. Scaling 
by 0.8 before foreground subtraction lowers the flux of the MIPS map by ^ 5000 Jy. 

The [CII] 158 /im line lies in the MIPS 160 /xm bandpass but not the nearby DIRBE bands. However, we show 
below that reasonable estimates of the contribution by this interstellar cooling line are too low to explain the excess 
flux in the MIPS map. Instead, the difference may be a result of the mapping method or simply an offset between the 
MIPS and DIRBE calibrations for the run over which the SMC data were taken. Throughout this paper we adopt the 
DIRBE absolute flux calibration, and use the rescaled MIPS map. We mention the effect that excess emission at 160 
/im would have on our results in the discussion of uncertainties. 

Figure [T4l shows the calibrated MIPS 160 /i,m data convolved to 0°.7 and 1'.5 and plotted against data from DIRBE 
and ISO. We bin the data by intensity and plot the mean and RMS scatter for each bin. Arbitrarily normalized 
histograms below each plot show the relative distribution of intensities for ISO and DIRBE and so give an idea of the 
importance for each point. The agreement between MIPS and the other data is good, particularly over the range of 
common intensities (high values in the histogram). The calibration could be improved by adding additional parameters 
(e.g. an intercept or quadratic term), but this runs the risk of skewing our results by overprocessing the data. Figure 
[13] may serve as a guide to the reader of the uncertainty in the calibration. 

In the remainder of this appendix, we show that the discrepancy between MIPS and DIRBE is too large to be 
accounted for by the [CII] 158 /im fine structure line or the known MIPS light leak. 

Is the Excess Due to [CII] Emission? 

The [CII] 158 /im line lies in the MIPS 160 /im bandpass. Does contamination by this important interstellar cooling 
line explain the discrepancy between the MIPS and DIRBE data? No extensive maps of the [CII] line in the SMC are 
available, so we estimate the [CII] contamination of the 160 /im map in two ways: by extrapolating from measurements 
o f [CII] in a few sta r -form ing r egions and by using the [CII]-to-FIR and [Cll]-to-CO ratio measured in other galaxies. 

[Israel fc MalonevI (|1993f l and lBolatto et all ()2000l) both measured t he [CI I] line towards several star-forming regions. 
Using the Kuiper Airborne Observatory (KAO). Israel fc MalonevI (|1993l ) found an integrated [Cll] line flux across 
five of the brightest star-forming regions that corresponds to ~ 50 Jy averaged across the 4 x 10^ Hz MIPS bandpass 
(F. Israel, private communication). Although the KAO measurements cover only a small portion of the SMC, they do 
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Fig. 14. — The MIPS 160 fim flux after correction (scaling by 0.8) and convolution to the DIRBE (left) and ISO (right) resolution plotted 
as a function of (left) DIRBE interpolated to 160 /im from nearby wavebands and (right) ISO. We bin the data by intensity and then plot 
the mean and RMS scatter in each bin. The agreement with ISO is good, but at high intensities the MIPS data is 20 - 30% lower than the 
DIRBE data, probably because of saturation effects. Normalized histograms show the distribution of intensities for both data sets. 



include the brightest star-forming regions which account for « 10% of the FIR emission. If [CII] emission is directly 
proportional to FIR emission across the SMC (rather than preferentially concentrated in PDR/star-forming regions) 
then [CII] emission from the SMC totals 500 Jy 

Staccy ct al. (1991) found that the [CII] emission accounts for 0.1 - 1 % of the FIR luminosity from a star-forming 
galaxy. The integrated FIR flux from the SMC based on IRAS i s 4 x 10~^° W m ~^ so we expect ^ 10^ - lO'^ Jy from 
the [CII] hue in the SMC based on the FIR emission. iMochizuki et al.l ()1998[ ) and lMadden et all (jlQQTl) report ratios 
between the [CII] and CO 1— >0 line fluxes of 10"* - 10^ for some actively star forming dwarf galaxies. iMizuno et al.l 
(|2001f ) finds the integrated CO intensity of the SMC to be ~ 10^ K km s~^ pc^, implying a flux from the CO line of 
3 X 10""'^'^ W m~^. If the SMC has a ratio of [CII] to CO emission consistent with the other dwarfs, this implies a [CII] 
luminosity of 10^ - 10^ Jy, consistent with the estimate from the FIR. 

Thus, several lines of argument suggest that the [CII] line contributes 10^ - 10"^ Jy to the MIPS 160 fim band. The 
most likely scenario is that a few hundred Jy of emission of total luminosity in the MIPS bandpass can be attributed 
to the [CII] line. This is significantly less than the ^ 5000 Jy needed to explain the discrepancy between the MIPS 
160 /xm band and the emission in nearby DIRBE bands. 

Is the Excess due to the MIPS Light Leak? 

The known light leak into MIPS from the near-IR cannot be the source of the flux discrepancy and does not represent 
a concern. According to the MIPS Data Handbook, data with a 160 /im / 2 /ini flux density ratio greater than 0.004 
will be uncorrupted by the light leak. Estimated from the DIRBE data, the median ratio of 160 /xm flux density to 2 
/im flux density along lines of sight towards the SMC is 75 and is never less than 35. 

PRODUCING THE DUST MAP 

In this section, we describe how we use the IRAS 100 /im and MIPS 160 /im data to construct a map of the dust 
mass surface density, ^dust, in the SMC. We then give brief descriptions of the sources of uncertainty that may affect 
this map. 

Method 

Deriving the dust mass surface density S^ust follows this basic method: 

1. Estimate the dust temperature, T, or distribution of temperatures the line of sight. 

2. Adopt an emissivity for the dust. This consists of a mass absorption coefficient, Ki/, and a power law wavelength 
dependence, cx A^'^, with an index /3. 

3. Convert the amplitude of the FIR emission into a dust mass surface density. For a single dust population: 



^Dust [g cm 
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[cm2 g-i] B^{T) [MJy ster"'] 
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where B^, (T) is the specific intensity of a blackbody with temperature T at a frequency v and is the intensity 
of measured FIR emission at the same frequency. 

Our data have hmited spatial resolution 75 pc at 100/im), and each line of sight through the SMC is probably 
several kiloparsecs long with more than one H I velocity component. As a result, we expected seve ral populat i ons o f 
dust to contri bute along each line of sight. Therefore, we treat the SMC using models presented bv lDale et al.l ()2001l ) 
and refined in lDale fc Heloul (|2002D . These models consider dust that is illuminated by a distribution of radiation fields, 
U, rang ing from 0.3 tolOf times the local intensity. They use an updated version of the description of interstellar 
dust by iDesert et al.l (|199C)D . The free parameters in the model are the power law index, a, of dust mass as a function 
of illuminating radiation field and the total dust mass, M^ust, proportional to the amplitude of the spectrum. These 
parameters, a and Mdust, are related by 



dMdust oc f/-" dU . (B2) 

So that the lower the value of a, the hotter the dust, on average. 

We derive T,dust for each line of sight by using the observed lOO-to-160 fim color to pick the best iDale fc Heloul 
(|2002) model, identified by its power law index a. From the amplitu de of the FIR emission and the model SED 
associated with that a we compute the dust mass along the line of sight. iDale et a l. (2001) emphasize that the models 
are well distinguished by the /eo/Aoo color alone. The lOO-to-160/xm color appears to do almost as good a job while 
avoiding the 24 - 70 fim bands. These bands introduce additional free parameters because they contain substa ntial 
contributions from polycyclic aromatic hydrocarbons (PAHs) and very small grains (VSGs) (|Desert et aLlfToOOh . the 
relative mix of which (i.e. the exact dust size distribution) is not well known in the SMC. 



Uncertainties 

Statistical: The statistical uncertainty associated with the map is modest. We repeatedly add noise with a magnitude 
from Table[2]to the data and estimate the la uncertainty in 'Sdust to be 8% at 2'. 6 resolution and 4% at 4' resolution. 
The uncertainty in Mdust is thus dominated by the choice of method, model, and systematic uncertainties in the data. 

Foreground Subtraction: The uncertainty associated with the foreground subtraction can have a sizeable impact, 
particularly in the diffuse (low intensity) Wing. We test the effect of errors in the foreground subtraction by adding 
and then subtracting 1 MJy ster~^ to the 100 fim map then the 160 /im map. This offset results in a ±20% change 
in the total dust mass over the SMC. In regions of low intensity an offset of 1 MJy ster^^ in one band may affect the 
dust surface density by up to 50%. 

Calibration and Saturation: If the MIPS calibration (rather than DIRBE) is correct then our MIPS 160 fim map 
should be scaled by ~ 1.25. This results in a dust content ~ 80% higher than we find. The possible saturation concerns 
for regions with /igo > 50 MJy ster~^ only have a small impact (< 10 %) on the deriv ed dust mass. 

Power Law vs. Single Population: The dust mass derived using the IDale fc Heloul (2002) models is a factor of ~ 3 
higher than one would derive using a single (3 = 2 population and the 100 and 160 /im bands to model each line of 
sight. Figure ITSl shows this difference by plotting the key quantity, T^dust per unit 160 /xm emission, as a function of 
the lOO-to-160 /xm color. Figure ITSl also shows the relative distribution of sight-lines with each color in the SMC and 
illustrates several other systematic effects discussed in this section. 

Umin cind Um.nx- The uppc r and l owcr limits to t he ra diation field distribution, Umin and Umax, may affect the 
derived dust mass. iDale et al] (|2001h : [Dale fc Heloul (|2002f) constructed their models to match the integrated SEDs of 
a sample of several dozen star-forming galaxies, so Umin = 0.3 and Umax — 10^ times the local value are motivated by 
observations. However, we test the effects of varying Umin and Umax on a simple model that includes only big grains 
with K250 — 8.5 cm^ g"-'^. If Umin = 0.1 is adopted instead of 0.3 then Mdust will be ~ 2 times higher. If Umin = 0.05 
instead of 0.3 then Mdust will be a factor of 4 higher. Higher values of Umin imply less dust mass, but Umin can not 
be much higher than « 0.5 or the models will not be able to reproduce the lowest /looZ-^ieo ratios seen in the SMC. 
Varying Umax has minimal effect over the range of colors observed in the SMC. Dust illuminated hy U = 0.3 has 
T « 14.5 K, so the uncertainty regarding Umin is a restatement of the potential population of very cold dust. Dust 
with T < 15 K has logiolioo/lieo < —0.5. Figure ITSi shows that there are very few lines of sight in the SMC with 
'o.gio/ioo 7-^160 < -0.3, let along -0.5. 

Choice of Bands (Inclusion of VSGs): The iDale fc Helo"ul (|2002f ) models account for the contribution from VSGs, 
albeit with a Galactic mix of grain types. Therefore, the choice to omit the 60 /im and 70 fim bands has a relatively 
small impact. Using the 100 and 160 fim data yields a mass 50 % higher than the 60 and 100 fim data, presumably 
because the SMC has more VSGs than the Milky Way mix adopted for the model. For the case of a single dust 
population it is critical to account for emission from VSGs if shorter wavelengths are included; failing to do so may 
result in estimates of dust mass es an order of magnitude too low (this ex plains a factor of ~ 5 out of the difference 
between our results and those of ISchwering|ll988t IStanimirovic et al.ll2000L the choice of a power law as opposed to a 
single model explains the balance). 

Mass Absorption Coefficient (n^): lAlton et al.l ()2004l ) synthesize determinations of the emissivity from comparisons 
of optical and sub-mm/FIR emission along the same line o f sight. They show that determinations of K250 span 
fror n fti 5 cm^ g~^ to 1 5 cm^ g~^, a substantial scatter. The 'Dale fc Hclou' (2002) models use a modifie d version of 
the iDesert et alj ()1990l ) dust model, which in turn adopts a slightly modified version of the iDraine fc Lei ri984) long 



wavelength mass absorption coefficient. Therefore the dust map in this paper adopts K250 ~ 8.5 cm^ g ^ (jDraine fc Le 
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TABLE 7 

Uncertainties in the Dust Map 
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^ Per line of sight. Negligible over the whole SMC.'' If the 
60 and 100 /im bands are used instead of the 100 and 160 
pm bands. In the limit, this is the same as the cold dust 
concern. 




Fig. 15. — A detailed look at the methodology used to create our dust map. (left) A normalized histogram of the lOO-to-160 /jm color 
for the SMC along with the dust mass surface density per unit 160 lira intensity, "i^duat/ Ii60i ^ ^ function of FIR color, logjQ /lOo/AeOi 
for several cas es. The t hree bold lines show the relationship for a single population of dust with f3 = 1.5, one with /3 = 2, and the power 
law models bv lDale fc H elou (2002). The thin lines show the same relationship for simple power law models with the minim u m (so lid) and 
maximum (dashed) radiation field varied. These lines show the impact of the assumptions that go into the IDale fc Heloul II2002I) model. 
The maximum radiation field has very little impact, but the minimum radiation field affects the conversion between FIR intensity and dust 
mass surface density, (right) A map of the FIR color, logj^Q /ioo/^i60i in the SMC. The star-forming regions have high lOO-to-160/im ratios 
and are easily identified from the map. 

I19841) . a widely used value close to the middle of the range of literature determinations. The dust surface densities we 
derive are uncertain by a factor of 2 due to the uncer tainty in the mass absorption coefficient. All of these values are 
not equally likely, but the value of K250 ~ 5 cm^ g"-'^ ()Li fc Drainel l200l'^ is also commonly used. 

Emissivity Dependence on Wavelength ((5): At long wavelengths, the wavelength dependence of Ky is usually treated 
as a power la w with an index /3 , so that (x A~ ^. However, ther e is substantial evidence that f3 may vary with 
environment ('Pup ae et all 120031 : [Alton et al.ll2004f) . The models by iDale fc Heloul (|2002D incorporate an em pirical 
radiation field de pend e nce fo r /3 calibrated using sub- mm data for the sample of galaxies used to develop the Dale et al.l 
(|200ll) : iDale fc'Heloul (j2002[ ) models. For this calibration f3 = 2.5, 2.1, and 1.3 for radiations fields 1, 10, and 1000 
times the local value. For the values of a we find in the SMC, /3 « 2. This may be too high, f3 = 1.5 may be more 
appropriate (see ^3.3\i . Because the 100 and 160 ^m bands are fairly close in wavelength, the impact of changing f3 by 
this amount is small; we would find Mdust to be « 20% lower for /3 = 1.5 than (3 = 2. Though it ha s a small effect on 
Mdust, P is important to interpreting the sub-mm excess observed towards the SMC (jAguirre et a l. 2003). Further, 
since P seems to be a function of environment, variations between diffuse and dense gas may systematically affect our 
H2 map (see main text). 



